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GREEN’S FUNCTION ASYMPTOTICS NEAR THE INTERNAL 
EDGES OF SPECTRA OF PERIODIC ELLIPTIC OPERATORS. 
SPECTRAL GAP INTERIOR. 

MINH KHA*, PETER KUCHMENT*^ AND ANDREW RAICH** 


Abstract. Precise asymptotics known for the Green function of the Laplacian 
have found their analogs for bounded below periodic elliptic operators of the second- 
order below and at the bottom of the spectrum. Due to the band-gap structure 
of the spectra of such operators, the question arises whether similar results can 
be obtained near or at the edges of spectral gaps. In a previous work, two of the 
authors considered the case of a spectral edge. The main result of this article is 
finding such asymptotics near a gap edge, for “generic” periodic elliptic operators 
of second-order with real coefficients in dimension d > 2, when the gap edge occurs 
at a symmetry point of the Brillouin zone. 


1. Introduction 

The behavior at inhnity of the Green function of the Laplacian in M” outside and 
at the boundary of its spectrum is well known. Analogous results below and at the 
lower boundary of the spectrum have been established for bounded below periodic 
elliptic operators of the second order in [3l|2l] (see also [3l] for the discrete version). 
Due to the band-gap structure of the spectra of such periodic operators, the question 
arises whether similar results can be obtained at or near the edges of spectral gaps. 
The corresponding result at the internal edges of the spectrum was established in [22] . 
The main result of this article, Theorem 12.111 is the description of such asymptotics 
near the spectral edge for generic periodic elliptic operators of second-order with real 
coefficients in dimension d > 2, if the spectral edge is attained at a symmetry point 
of the Brillouin zone. 

It is well known that outside of the spectrum the Green function decays exponen¬ 
tially at inhnity, with the rate of decay controlled by the distance to the spectrum. 
See, e.g., well known Gombes-Thomas estimates laizi- However, comparison with 
the formulas for the case of the Laplacian shows that an additional algebraically de¬ 
caying factor (depending on the dimension) is lost in this approach. Moreover, the 
exponential decay in general is expected to be anisotropic, while the operator theory 
approach can provide only isotropic estimates. The result of this paper provides the 
exact principal term of asymptotics, thus resolving these issues. 
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2. Assumptions, notation and the main result 
C onsider a linear second order elliptic operator in with periodic coefficients 

d 

(1) L{x, D)=Y, Dk{au{x)Di) + V{x) = D*A{x)D + V{x). 

k,l=l 

d 

Here A = (aki)i<k,i<d, D = (Di,..., and Dk := -idk = coef¬ 

ficients ttkh V are smooth real-valned fnnctions on periodic with respect to the 
integer lattice in i.e., au{x+n) = (x) and V(x+n) = V(x),Vx G n G 

The operator L is assnmed to be elliptic, i.e., the matrix A is symmetric and 

d 

(2) aki{x)^k^i > 

k,l=l 

for some 9 > 0 and any x G ^ The operator L, with 

the Sobolev space as the domain, is an nnbonnded, self-adjoint operator in 

(see e.g., [29]). 

The spectrnm of the above operator L in has a band-gap structure 

imgEB], i.e., it is the nnion of a sequence of closed bounded intervals (bands or 
stability zones of the operator L) [aj, Pj] C M (j = 1, 2,...): 

OO 

(3) a{L) = [j[aj,Pj], 

i=i 

such that aj < cij+i, (3j < (3j+i and limj^oo ctj = oo. The bands can (and do) overlap 
when d > 1, but they may leave open intervals in between, called spectral gaps. 
Thus, a spectral gap is an interval of the form for some j G N for which 

Q!j+i > f3j. We make a convention that the open interval (— cxo, oi), which contains all 
real numbers below the bottom of the spectrum of L, is also a spectral gap. However, 
we will be mostly interested in finite spectral gaps. 

In this text, we study Green’s function asymptotics for the operator L in a spectral 
gap, near to a spectral gap edge. More precisely, consider a finite spectral gap 
{/3j,aj+i) for some j G M and a value A G {(3j,aj+i) which is close either to the 
spectral edge f3j or to the spectral edge Uj+i. We would like to study the asymptotic 
behavior when |x — |/| —)■ cxo of the Green’s kernel G\{x,y) of the resolvent operator 
Rx^l ■= {L — A)“^. The case of the spectral edges (i.e., A = cij+i or A = /3j) 
was studied for the similar purpose in [22]. All asymptotics here and also in [22] 
are deduced from an assumed “generic” spectral edge behavior of the dispersion 
relation of the operator L, which we will briefly review below. 

Let W = [0,1]'^ C be the unit cube, which is a fundamental domain of 
with respect to the lattice Z'’* (Wigner-Seitz cell). The dual (or reciprocal) 
lattice is 27rZ'^ and its fundamental domain is [—7r,7r]'^ (Brillouin zone). 

The d-dimensional tori with respect to the lattices Z'’* and 27rZ'’* are denoted by 
T := and T* := M'^/27rZ'’*, respectively. 
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Definition 2.1. For any k G C'^, the subspace H^iW) C consists of restric¬ 
tions to W of functions / G that satisfy for any 7 G the Floquet-Bloch 

condition (also known as automorphicity condition or cyclic condition) 

(4) /(x + 7 ) = f{x) for a.e x G IF. 

Here denotes the standard Sobolev space of order s. Note that when s = 0, the 
above space coincides with LP'iW) for any k. 

Due to periodicity, the operator L(x, D) preserves condition (|1|) and thus, it defines 
an operator L{k) in L^(IF) with the domain HliW). In this model, L{k) is realized 
as a /c-independent differential expression L{x, D) acting on functions in W with 
boundary conditions depending on k (which can be identified with sections of a 
linear bundle over the torus T). An alternative definition of L{k) is as the operator 
L(x, D + k) in L^(T) with the domain In the latter model, L{k) acts on the 

fc-independent domain of periodic functions on W as follows: 

(5) e-*^'^L(x, = {D + kyA{x){D + k) + V{x). 

We use the latter model of L{k) throughout this paper, unless specified differently. 

Note that the condition (0]) is invariant under translations of k by elements of the 
dual lattice 27rZ'^. Moreover, the operator L{k) is unitarily equivalent to L(fc-|-27r7), 
for any 7 G Z'^. In particular, when dealing with real values of k, it suffices to 
restrict k to the Brillouin zone [—vr, tt]'^ (or any its shifted copy). It is well-known 
(see [SIDS1I2S] ) that the spectrum of L is the union of all the spectra of L{k) when k 
runs over the Brillouin zone, i.e. 

(6) a{L) = IJ (T{L{k)). 

k£[—7T,7r]^ 

Hence, the spectrum of L is the range of the multivalued function 

(7) k ^ X{k) := a{L{k)), fc G [—vr, vr]'^. 

By ([2]), L{k) is self-adjoint in L^(T) and has domain i7^(T) for each k G By 
the ellipticity of L, each L{k) is bounded from below and has compact resolvent. This 
forces each of the operators T(fc), k G to have discrete spectrum in M. Therefore, 
we can label its eigenvalues in non-decreasing order: 

( 8 ) Xi{k) < X2{k) < ... . 

Hence, we can single out continuous and piecewise-analytic band functions Xj{k) 
for each j G N [53] • The range of the band function Xj constitutes exactly the band 
[aj^(3j\ of the spectrum of L shown in (|5|). 

Definition 2.2. A Bloch solution of the equation L(x, D)u = 0 is a solution of the 
form 

n(x) = e^’^'^(j){x), 

where the function 0 is 1-periodic in each variable xj for j = 1,... ,d. The vector k 
is the quasimomentum and z = = (e**'i,..., is the Floquet exponent (or 

Floquet multiplier) of the solution. In our formulation, allowing quasimomenta k 
to be complex is essential. 
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Definition 2.3. The (complex) Bloch variety of the operator L consists of all 
pairs (fc, A) G such that the equation Lu = Xu in has a non-zero Bloch 

solution u with a quasimomentum k. Similarly, the real Bloch variety i?L,R is 
Bl nM'^+h 

The Bloch variety can be treated as the dispersion relation/curve, i.e., the 
graph of the multivalued function \{k): 

BL = {{k,X):Xea{L{k))}. 

Note that L{k) is non-self-adjoint if A; ^ However, L{k) — L{0) is an operator of 
lower order for each k G C'^. Therefore, the spectra of all operators L{k) on the torus 
T are discrete (see pp. 188-190 in [T]). 

Remark 2.4. In fact, the main techniques of Floquet theory (e.g., ([6])) apply to non¬ 
self-adjoint operators. It is required only that the operators L{k) = L(a;, D + k) ■. 

—>■ L^(T) are Fredholm for k G C'^. The latter condition is always satisfied 
due to ellipticity and embedding theorems (see Theorem 2.1 in [H]). Unlike the 
self-adjoint case though, we do not have the band-gap structure as in (|3]). 

Definition 2.5. The (complex) Fermi surface F^ x of the operator L at the energy 
level A G C consists of all quasimomenta fc G C'^ such that the equation Lu = Xu in 
has a non-zero Bloch solution u with a quasimomentum k. For A = 0, we simply 
write Fl instead of Fi^. The real Fermi surface is Fl n 

Equivalently, k G Fl^\ means the existence of a nonzero periodic solution u of the 
equation L{k)u = Xu. In other words, Fermi surfaces are level sets of the dispersion 
relation. 

The following result can be found in Theorem 3.1.7 in [T9] : 

Lemma 2.6. There exist entire -periodic in k) functions of finite orders on 
and on such that the Fermi and Bloch varieties are the sets of all zeros of these 
functions respectively. 

From this lemma and the proof of Lemma 4.5.1 in [19] (see also [33]), the band 
functions X{k) are piecewise analytic on C*^. 

From now on, we fix L as a self-adjoint elliptic operator of the form ([1]), whose 
band-gap structure is as ([3|). By adding a constant to the operator L if necessary, we 
can assume that the spectral edge of interest is 0. It is also enough to suppose that 
the adjacent spectral band is of the form [0, a] for some a > 0 since the case when 
the spectral edge 0 is the maximum of its adjacent spectral band is treated similarly. 

Suppose there is no spectrum for small negative values of A and hence there is a 
spectral gap below 0. Thus, there exists at least one band function Xj{k) for some 
j G N such that 0 is the minimal value of this function on the Brillouin zone. 

To establish our main result, we need to impose the following analytic assumption 
on the dispersion curve Xj as in [22] : 


Assumption A 
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There exists G [—7r,7r]'^ and a band function Xj{k) such that: 

A1 Xj{ko) = 0. 

A2 \Xiik)\ > 0. 

A3 ko is the onl'^ (modulo 27rZ'^j minimum of Xj. 

A4 Xj{k) is a Morse function near ko, i.e., its Hessian matrix H := Hess {Xj){ko) 
at ko is positive definite. In particular, the Taylor expansion of Xj at ko is: 

= ^(^ - kofH{k - ko) + 0 {\k - ko\^). 

It is known la that the conditions A1 and A2 ‘generically’ hold (i.e., they can 
be achieved by small pertnrbation of coefficients of the operator) for Schrodinger 
operators. Although this has not been proven, the conditions A3 and A4 are widely 
believed (both in the mathematics and physics literature) to hold ‘generically’. In 
other words, it is conjectured that for a ‘generic’ selfadjoint second-order elliptic op¬ 
erator with periodic coefficients on each of the spectral gap’s endpoints is a unique 
(modulo the dual lattice nondegenerate extremum of a single band function 

Xj{k) (see e.g.. Conjecture 5.1 in [21] )• It is known that for a non-magnetic pe¬ 
riodic Schrodinger operator, the bottom of the spectrum always corresponds to a 
non-degenerate minimum of Ai |T 6 |. A similar statement is correct for a wider class 
of ‘factorable’ operators W- The following condition on ko will also be needed: 

A5 The quasimomentum ko is a high symmetry point of the Brillouin zone, 

i.e., all components of ko must he either equal to 0 or to n. 

We denote by X the set of such high symmetry points in the Brillouin zone. 

It is known [T3| that the condition A5 is not always satished and spectral edges 
could occur deeply inside the Brillouin zone. However, as it is discussed in [Tl|, in 
many practical cases (e.g., in the media close to homogeneous) this condition holds. 

We would like to introduce a suitable fundamental domain with respect to the dual 
lattice 27rZ'‘* to work with. 

Definition 2.7. Consider the quasimomentum ko in our assumptions. By A5, 
ko = (( 5 i 7 r, 6271 ,..., SdTi), where 6 j G {0, 1} for j G {!,..., d}. We denote by O 
the fundamental domain so that ko is its center of symmetry, i.e., 

d 

i=i 

When ko = t), O is just the Brillouin zone. 

We now introduce notation that will be used throughout the paper. 

^Finitely many such points can be also easily handled. 
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Notation 2.8. 

(a) Let Zi G C, 2^2 e z^, G and r* be positive numbers for i = 1,2,3. Then 

we denote by i?(zi,ri), D'{z 2 ,r 2 ) and D{z^,r^) the open balls (or discs) centered 
at zi, Z 2 and Z 3 whose radii are ri, r 2 and rs in C, and respectively. 

(b) The real parts of a complex vector z, or of a complex matrix A are denoted by 
3 ^( 2 :) and 3fJ(A) respectively. 

(c) The standard notation 0{\x — y\~'^) for a function / defined on means there 
exist constants C > 0 and R > 0 such that \f{x,y)\ < C\x — y\~'^ whenever 
\x — y\ > R. Also, f{x,y) = o{\x — y\~'^) means that 

lim |/(a;,i/)|/|x-|/|" = 0. 

\x—y\^oo 

(d) We often use the notation A < to mean that the quantity A is less or equal 
than the quantity B up to some multiplicative constant factor, which does not 
affect the arguments. 

As we discussed, for each 2 ; G C'^, the operator L{z) has discrete spectrum and 
is therefore a closed operator with non-empty resolvent set. These operators have 
the same domain and for each (j) G iL^(T), L{z)(j) is a L^(T)-valued analytic 

function of z, due to ([5]). Consequently, {L{z)}z^cd' is an analytic family of type A 
in the sense of Kato pisj^ . Due to A1-A2, Xj{ko) is a simple eigenvalue of L{ko). 
By using analytic perturbation theory for the family {L{z)}^^c<^ (see e.g.. Theorem 
XII.8 in [2H]), there is an open neighborhood V of k^ in and some Cq > 0 such 
that 

(PI) Xj is analytic in a neighborhood of the closure of V. 

(P2) \j(z) has algebraic multiplicity one, i.e., it is a simple eigenvalue of L{z) for 
any z G V. 

(P3) The only eigenvalue of L(z) contained in the closed disc 5(0, Cq) is Xj(z). 
Moreover, we may also assume that |Aj( 2 :)| < eo for each z eV. 

(P4) For each z E V, let (j){z, x) be a nonzero Z'^-periodic function of x such that 
it is the unique (up to a constant factor) eigenfunction of L{z) with the eigenvalue 
Aj(z), i.e., L{z)(l){z, •) = \j{z)(j){z, •). We will also use sometimes the notation 0^ for 
the eigenfunction 0(z, •). 

By elliptic regularity, (j){z,x) is smooth in x. On a neighborhood of V, (j){z, •) is a 
5^(T)-valued holomorphic function. 

(P5) By condition A4 and the continuity of Hess (Aj) we can assume that for all 
^ G H, 

23f?(Hess (Aj)( 2 ;)) > min (j(Hess (Aj)(/co))/dxd- 
(P6) V is invariant under complex conjugation. Furthermore, the smooth function 

(9) 5 ( 2 ;) := (0(z, •),0(T, •))l2(t) 

is non-zero on V, due to analyticity of the mapping 2 ; 1 —)■ (hlz, ■) and the inequality 

^It is also an analytic family in the Banach space of bounded linear operators acting from 77^ (T) 
to L2 (t). 

^The Hessian matrix of Xj. 
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The following lemma will be useful when dealing with operators having real and 
smooth coefficients: 

Lemma 2.9. (i) For k in and i G M, 

(10) K{k) = K{-k). 

In other words, each band Aj of L is an even function on 

(a) If ko G X, we have Xi{k + k^) = Xi{—k + ko) for all k in and i G N. 

Proof. Let <fk be an eigenfunction of L{k) corresponding to Xj{k). This means that 
(fk is a periodic solution to the equation 

( 11 ) L{x, d + ik)(j)k{x) = Xj{k)(j)k{x). 

Taking the complex conjugate of flTT]) . we get 

L{x, d - ik)(j)k{x) = Xj{k)(f)k{x). 

Therefore, (fk is an eigenfunction of L{—k) with eigenvalue Xj{k). This implies the 
identity (HU . 

(ii) By (i), Xi{k + ko) = Xi{—k — ko) = Xi{—k + ko) since 2ko G 27rZ‘^. □ 

Corollary 2.10. ///9 G such that kfj + i(5 E V then Xj{ko + ijl) G M. 

Proof. Indeed, the statement (ii) of Lemma [2.91 implies that the Taylor series of X{k) 
at ko has only even degree terms and real coefficients. □ 

Corollary I2.1UI allows us to define near /3 = 0 the real analytic function E{/3) : = 
Xj{ko + i/3) near 0. Since its Hessian at 0 is negative-definite (by A4), there exists 
a connected and bounded neighborhood Vq of 0 in such that ko + iVo C V and 
Hess {E){/3) is negative-definite whenever (3 belongs to Vq. Thus, E is strictly concave 
on Vo and sup^gy^ E{I3) = i?(0) = 0, VE{/3) = 0 iff = 0. Note that at the bottom 
of the spectrum (i.e., j = 1), we could take Vq as the whole Euclidean space 

By the Morse lemma and the fact that 0 is a nondegenerate critical point of E, there 
is a smooth change of coordinates </) : Uq ^ so that 0 E Uq CC Vq, Uq is connected, 
0(0) = 0 and E(0“^(a)) = —|ap,Va G 4>{Uo). Set Kx := {[3 E Uq '■ E{f3) > A} and 
Ta := {/d G Ho : E{/3) = A} for each A G M. Now, we consider A to be in the set 
{ —|ap : a E 0(t/o),a 7 ^ 0}. Then Kx is a strictly convex d-dimensional compact 
body in R'^, and Ta = dKx is a compact hypersurface in R*^. The compactness of Kx 
follows from the equation —|0(/3)p = E{/3) > X which yields that \/3\ = |0“^(0(/3))| < 
max{|0“^(a)| : a E 0(Ho), |ap < —A}. Additionally, limA^.o- max^g^^ |/3| = 0. 

Let /Ca be the Gauss-Kronecker curvature of Ta- Since the Hessian of E is negative- 
definite on Ta, /Ca is nowhere-zero. For the value of A described in the previous 
paragraph and each s E there is a unique vector G Ta such that the value of 
the Gauss map of the hypersurface Ta at this point coincides with s, i.e. 

(12) VE{l3s) = -\VE{l3s)\s. 

This is due to the fact that the Gauss map of a compact, connected oriented hyper¬ 
surface in R'^, whose Gauss-Kronecker curvature is nowhere zero, is a diffeomorphism 
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onto the sphere ^ (see e.g., Theorem 5, p.l04 in |32] or Corollary 3.1 in mi)- 
Thus, (3s depends smoothly on s. We also see that 

lim max |/3s| = 0. 

|A |^0 

Note that could be dehned equivalently by using the support functional h of the 
strictly convex set K\. Recall that for each direction s G 

h(s) = max(s, 

Then /3s is the unique point in Ta such that {s,(3s) = h{s). 

By letting |A| close enough to 0, we can make sure that (— G Then 

(13) {ko + it(3s, (t, s) G [0,1] X c V. 

We can now state the main result of the paper. 


Theorem 2.11. Suppose conditions A1-A5 are satisfied. For A < 0 sufficiently 
close to 0 (depending on the dispersion branch Xj and the operator L), the Green’s 
function G\ of L at X admits the following asymptotics as \x — y\ ^ oo: 


(14) 


Gx{x,y) = 


^{x-y)(iko-l3s) 


|V.B(/3,)|('^-3)/2 


0fco+j/3s {^)4’ko—i0s {y) 


{27i\x - y\)G-^)G fjggg {E){/3s)VsY^^ (0fco+i/3s, 0fco-i/3jihT) 


+ e 


{y-x).y. 


r{x,y). 


Here s = {x — y)/\x — y\, Vs is the projection from onto the tangent space of the 
unit sphere at the point s, and for any e > 0, there exists a constant > 0 
(independent of s) such that the remainder term r satisfies \r{x,y)\ < Cs\x — y\~‘^V+^ 
when \x — y\ is large enough. 


This result achieves our stated goal of showing the precise (anisotropic) rates of 
the exponential decay of the Green’s function and capturing the additional algebraic 
decay factor. 


3. Proof of the main theorem 12.111 and some remarks 


Theorem 12.111 is a direct consequence of its local (with respect to the direction of 
{x — y)) version: 


Theorem 3.1. Under the hypotheses of Theorem. \2.11\ and when A ~ 0, for each 
u G there exists a neighborhood in containing oj and a smooth function 
e(s) = {es^ 2 , ■ ■ ■, Gs,d) ■ Vij —)■ {TsE)'^~^Y~^ such that the asymptotics 


Gx{x,y) = 


(15) 


|VE(/?,)|(''-3) 


^(x-y){iko-ys) 

{2'k\x - y\)G-TG \^det (- 6 ^,^ • Hess {,E){j3s)es,q)2<p,q<d 

^y~’’’^'^‘’r{x,y), 


1/2 


4^ko+il3s{^)4^ko-iysiy) I „ 

X —;-;-^o 


ifiko+ilSs^ 4>ko-ilSs)L^{T) 


hold for all {x,y) such that s = {x — y)/\x — y\ G V^j. Furthermore, \r{x,y)\ < 
C{e,u)\x — y\~F‘^+^ for some positive constant C{e,u) depending on e and u. 
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Proof of Theorem 2.11. Observe that for any basis {esi} 2 <i<d of the tangent space 

det {-Vs Hess {E){l3s)Vs) = det (-e^^p • Hess {E){l3s)es,q)2<p,q<d- 

Now, using of a hnite cover of the unit sphere by neighborhoods Vu)^ in Theorem I3.H 
one obtains Theorem 12.111 □ 


Remark 3.2. 


• The asymptotics (ITTll (or (|T5ll ) resemble the formula (1.1) in Theorem 1.1 [2l] 
when A is below the bottom of the spectrum of the operator. Moreover, as in 
Theorem 1.1 in [25], by using the Gauss-Kronecker curvature /Ca, the main 
result (fTT|) could be restated as follows: 


^ _^_ (l)kq+ipXx)(fkq-ifiXy) 

(2vr|a:-2/|)G-b/2 \VE{^)\lC^{^yV (0^0+*/^., 

yy-vp.o{\x - 


• Although flTTl) is an anisotropic formula, it is not hard to obtain from flTTll an 
isotropic upper estimate for the Green’s function G\ based on the distance 
from A to the spectrum of the operator hj^e.g., there are some positive 
constants Gi,G 2 (depending only on L and \j) and (which may depend 
on A) such that whenever \x — y\ > G 3 , the following inequality holds: 




-C2|A|1/2|3;_p| 


\x — 

If the band edge occurs at hnitely many points, rather than a single ko, one 
just needs to combine the asymptotics coming from all these isolated minima. 


Now we outline the proof of Theorem 13.11 In Section 5, we introduce the tools 
of Floquet-Bloch theory to reduce the problem to that of hnding the asymptotics 
of a scalar integral, similarly to [22]. The purpose of Section 4 is to prepare for 
Section 5, i.e., to shift the integral from the fundamental domain O along some 
purely imaginary directions in C'^. Section 6 is devoted to estimating the leading 
term integral by adapting the method similar to the one used in the discrete case 
[34] . The decay of the remainder r{x,y) comes from some elementary estimates of 
the difference between the scalar integral and the leading term of its asymptotics. In 
order to not overload the main text with technicalities, the proofs of some auxiliary 
statements are postponed till Sections 7-9. 


4. On local geometry oe the resolvent set 

The following proposition shows that for any s G fco + ^/ds is the only complex 

quasimomentum having the form of A: -|- it^s where k ^ 0,t ^ [0,1] such that A is 
in the spectrum of the corresponding hber operator L{k + it/3s). In other words, by 
moving from fc G O in the direction iflg, the hrst time we hit the Fermi surface 
El^x (i.e., the spectrum of L{k) meets A) is at the value of the quasimomentum 

^Recall that the spectral edge is assumed to be zero. 
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k = ko + ij3s- This step is crucial for setting up the scalar integral in the next section, 
which is solely responsible for the main term asymptotics of our Green’s function. 

Proposition 4.1. //|A| is small enough (depending on the dispersion branch Xj and 
L), then A G p{L{k + itfds)) if and only if {k, t) ^ (fco, !)■ 

The proof of this proposition is presented in Subsection 19.31 


5. A Floquet reduction of the problem 


We recall here some basic properties of the Floquet transform and then apply 
this transform to reduce our problem to finding asymptotics of a scalar integral 
expression, which is close to the one arising when dealing with the Green’s function 
of the Laplacian at a small negative level A. As in [22], the idea is to show that only 
the branch of the dispersion relation Aj appearing in the Assumption A dominates 
the asymptotics. 


5.1. The Floquet transform. Let us consider a sufficiently fast decaying function 
f{x) (to begin with, compactly supported functions) on W^. We need the following 
transform that plays the role of the Fourier transform for the periodic case [T^[26]. 
In fact, it is a version of the Fourier transform on the group Z'’* of periods. We use 
the following version, which is slightly different from the one used in [22] . 


Definition 5.1. The Floquet transform T 


f{x) f{k,x) 


maps a function / on into a function / defined on x Mf' in the following way: 


f{k,x) := + 

'yeZ'^ 


From the above definition, one can see that / is ^Aperiodic in the x-variable and 
satisfies a cyclic condition with respect to k: 

i f{k,x + j) = f{k,x), Vy e 
\ f{k + 27r7, x) = x), Vy G Z'^ ' 

Thus, it suffices to consider the Floquet transform / as a function defined on (9 x T. 
Usually, we will regard / as a function f{k, •) in fc-variable in O with values in the 
function space L^(T). 

For our purpose, we need to list some well-known results of the Floquet transform 
(see e.g., [19]): 

Lemma 5.2. 1. The transform is an isometry of onto 








and of into 
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2. The inversion T ^ is given by the formula 

(16) f{x) = (271)-'^ [ e^^-^f{k,x)dk, X e W^. 

Jo 

By using cyclic conditions of f, we obtain an alternative inversion formula 

(17) f{x) = {2'K)-<^ ! e^^-^f{k,x--i)dk, xeiy + 7. 

Jo 

3. The action of any TJ-periodic elliptic operator L (not necessarily self-adjoint) 
in L^(R‘^) under the Floguet transform 3F is given by 

r® r® 

J^L{x,D)J^^^ = / L{x,D + k)dk= / L{k)dk, 

Jo Jo 

where L{k) is defined in 
Eguivalently, 

rf{k) = mm, v/ e h^r^). 

4 . (A Paley-Wiener theorem for P.) Let (j){k,x) be a function defined on 
X R^ such that for each k, it belongs to the Sobolev space H^{T) for s G R"*" and 

satisfies the cyclic condition in k-variable. Then 

(1) Suppose the mapping k —)■ 0(/c, •) is a C°°-map from R'^ into the Hilbert space 
7/^(T). Then (j){k,x) is the Floguet transform of a function f G H^{R'^) such 
that for any compact set K in R'^ and any N > 0, the norm \\f\\Ho{K+'y) < 
Cn\' 1\~^■ In particular, by Sobolev’s embedding theorem, if s > d/2, then the 
pointwise estimation holds: 

|/(t)| <C'^(1 + |t|)-'^, ViV>0. 

(2) Suppose the mapping k (j){k, •) is an analytic map from R^ into the Hilbert 

space i7®(T). Then (j){k, x) is the Floguet transform of a function f G i7^(R'^) 
such that for any compact set K in R'^, one has \\f\\H‘{K+'y) < In 

particular, by Sobolev’s embedding theorem, if s > d/2, then the pointwise 
estimation holds: 

5.2. The Floquet reduction. 

The Green’s function G\ of L at A is the Schwartz kernel of the resolvent operator 
R\ = {L — A)“^. Fix a A < 0 such that the statement of Proposition 14.11 holds. For 
any s G and t G [0,1], we consider the following operator with real coefficients 
on R'’*: 

(18) Li,, := 

For simplicity, we write L, := Li,, and note that Lq,, = L. Due to self-adjointness 
of L, the adjoint of Li,, is 

( 19 ) LI, = 
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By definition, Lt^sik) = L{k + it(3s) for any k in C'^ and therefore, ([6]) yields 

(20) cr{Lt^s) = IJ a{L{k + it^s)) ^ {Xj{k + itl3s)}k&o- 

keo 

The Schwartz kernel of the resolvent operator Rg x := (L^ — A)“^ is 

(21) G,,A(a:, y) = e^-"GA(x, i/)e-^-^ = e^^X^-y^Gx{x, y). 

Thus, instead of hnding asymptotics of G\, we can focus on the asymptotics of Gg^x- 
By fl2ni) and Proposition 14.11 A is not in the spectrum of Lt^g for any s G and 
t G [0,1). Let us consider 

R,.s.xf - (iM - W‘/. / e iL„p(K''). 

where stands for compactly supported functions in 

Applying Lemma [5.21 we have 

R^^fik) = {Lt,g{k) - X)-^f{k), (f, k) G [0,1) X (P. 

We consider the sesquilinear form 

(^M,a/,^) = {{LUk) - X)-^m,m) dk, 

where cp G 

In the next lemma (see Subsection 9.2 for its proof), we show the weak convergence 
of Rt^sp in L^omp as f 1 and introduce the limit operator Rg x = lim Rt^gp. The 

t —^ 1 

limit operator Rg^x is central in our study of the asymptotics of the Green’s function. 
Lemma 5.3. Let d > 2. Under Assumption A, the following equality holds: 

( 22 ) \^_{Rt,g,xf,g:>) = ( 27 r)-"^ (^.(fc) - X)-^f{k),tp{k)) dk. 

The integral in the right hand side of absolutely convergent for f, ip in L^omp(®'^)- 

Thus, the Green’s function Gg^x is the integral kernel of the operator Rg^x defined as 
follows 

(23) ^f{k) = iLgik)-X)-^fik). 

5.3. Singling out the principal term in Rgp. By fl23H . the Green’s function Ggp 
is the integral kernel of the operator Rg^x with the domain L^^^p(M'’*). The inversion 
formula flT^ gives 

Rg,xf{x) = ( 271 )-" [ U’^-^{Lg{k) - X)-^f{k, x)dk, X G M". 

Jo 

The purpose of this part is to single out the part of the above integral that is respon¬ 
sible for the leading term of the Green’s function asymptotics. 

To hud the Schwartz kernel of Rgp, it suffices to consider functions / G G^(M"). 
Our hrst step is to localize the integral around the point fco- Let us consider a 
connected neighborhood V of ko on which there exist nonzero Z"-periodic (in x) 
functions (j)z{x),z G V satisfying 1) L{z)(j)z = Xj{z)(f)z and 2) each (fz spans the 
eigenspace corresponding to the eigenvalue Xj{z) of the operator L{z). According to 
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(P3), \j{y) C 5(0, eo) and 95(0, eo) C p[L{z)) when z &V. For snch z, let P{z) 
be the Riesz projection of L{z) that projects L^(T) onto the eigenspace spanned by 
(pzi i-G., 

P{z) = —TT— (f {L{^) ~ oi)~^da. 

J |a|=eo 

Taking the adjoint, we get 

P{z)* = i {L{z) - a)-^da = P{z), 

which is the Riesz projection from L^(T) onto the eigenspace spanned by (f)^- Recall 
that dne to o, by choosing |A| small enongh, there exists ro > 0 (independent of 
s) snch that fc ± i/d* G R for fc G D{ko,ro) fl We denote Ps{k) := P{k + for 
snch real k. Then Ps{k) is the projector onto the eigenspace spanned by (j){k + i[ds) 
and Ps{k)* = P{k — i/3s). Additionally, dne to (P6), 


(24) 


Ps{k)g 


{9,4>{k-i/3s))LHT} 

{cj){k + i/3s), (f){k - */3s))l2(t) 




Wg G L\T). 


Let ?7 be a cnt-off smooth fnnction on O snch that snpp(? 7 ) d D{ko,ro) and rj = 1 
aronnd k^. 

We decompose / = ? 7 / + (l — 77 )/. When k 7 ^ k^, the operator Ls{k) — \ is invertible 
by Proposition 14.11 Hence, the following fnnction is well-defined and smooth with 
respect to (/c, x) on x 


Ug{k, x) = {Ls{k) - A) ^(1 - g{k))f{k, x). 


Using Lemma 15.21 smoothness of Ug implies that Ug has rapid decay in x. Now we 
want to solve 


(25) {Ls{k) - X)u{k) = g{k)f{k). 

Let Qa{k) = I — Ps{k) and we denote the ranges of projectors Ps{k), Qs{k) by 
R{Ps{k)), R{Qs{k)) respectively. We are interested in decomposing the solution u into 
a sum of the form U 1 +U 2 where Ui = Ps{k)ui and U 2 = Qs{k)u 2 - Let fi = Ps{k)g{k) f 
and /2 = Qs{k)g{k)f. Observe that since the Riesz projection Ps{k) commutes with 
the operator Ls{k) and R{Ps{k)) is invariant under the action of Ls{k), we have 
Qs{k)Ls{k)Ps{k) = Ps{k)Ls{k)Qs{k) = 0 and Qs{k)Ls{k)Qs{k) = Ls{k)Qs{k). Thus, 
the problem of solving fl25D can be reduced to the following block-matrix structure 
form 


{La{k) - X)Ps{k) 

0 

0 

{La{k) - X)Qa{k) 


Ml 

7*2 



When k is close to ko, 


5(0, eo) n a(L^(A;)|ij(Q^(fc))) = 5(0, Cq) O a{L{k + */?J) \ {\j{k + = 0. 

Since A = Aj(fco + */ds) G 5(0, cq), A must belong to p(Ls(fc)|R(Q 4 fc)))- Hence, the 
operator function U2{k) = {Ls{k) — X)~^Qs{k)f2{k) is well-dehned and smooth in k 
and hence by Lemma [5.21 again. U 2 has rapid decay when |x| —)■ 00 . Indeed, we claim 
that the Schwartz kernel coming from the operator-valued function {l — r]{k)){Ls{k) — 
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X)~^ + ri{k){{Ls{k) — X)\ji(^Q^(^k)))~^Qs{.k) decays fast enough to be included in the error 
term r{x,y) in ffl^ . We shall give a microlocal proof of this claim in Section 7. 

The Ml term contributes the leading asymptotics for the Schwartz kernel Gs,x- 
Therefore, we only need to solve the equation {Ls{k) — X)Pa{k)ui = fi on the one¬ 
dimensional range of Ps{k). 

Applying fl2T|) . we can rewrite 


m 


r]{k){f,(j){k-i/3s))LHT) ,,, .. . 


so that equation becomes 


{La{k)-X) 


(Mi,0(fc-Z/3.)U2(T) 


V{k){l<t^{k-Z^a))LHT) 


So, 


{Xj{k + i^a) - A)(Mi, 0(fc - iPa))LHT) , i h(^)(/, “ */^s))l2(T) - o , ^ 

(cjz{k + zfia),(^{k-Z/3a))LHT) ^ + Z^a), ^{k - Z/3a)) L^T)^ 

In addition to the equation ui = Pa{k)ui, Ui must also satisfy 


{Xj{k + if3a) - X){ui, (j){k - if3a))L^(T) = r]{k){f, (t){k - i/3J)L2(T)- 
Thus, we define 

r?(fc)0(fc -h Zl3a, •)(/, 0(fc - i/3s))L2(T) 

’ ’ {(t){k + zl3a),(t){k-zl3a))L^(T){Xj{k + zl3a)-Xy 

By the inverse Floquet transform fITBD . 


Mi(a;) = ( 27 r) 


-d f ^ik-x V{k)(j){k + i/3a,x){f,y{k - i/3a))L^T) 

Jo^ {(t){k + i(3a),(t){k - z(3a))L^{j){Xj{k P i(3a) - A) 


dk, 


for any x G 


5.4. A reduced Green’s function. We are now ready for setting up a reduced 
Green’s function Go, whose asymptotic behavior reflects exactly the leading term 
of the asymptotics of the Green’s function Ga,\. We introduce Go{x,y) (roughly 
speaking) as the Schwartz kernel of the restriction of the operator Rg^x onto the 
one-dimensional range of Pg (which is the direct integral of idempotents Pa{k)) as 
follows: 

Ui{x) = [ Go{x,y)f{y)dy, x 
Jm.d' 


where / is in L2^^p(M'^). 
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We recall from ([9]) that F{k + i(3s) is the inner product {(j){k + i(5a),4’{k — i(3s))L^{T)- 
As in [22], we notice that 


ui{x) = (27r) 


-d 





O JT 


F{k + i(5a){Xj{k + il3s) - A) 


= (2vr) / r]{k) 


[ 0 , 1 ]'' 


= (2vr)-'^ 


'O 


7eZ'' 


/[o,ip+7 F(fc + z/?,)(A,(A; + ^/3,)-A) ^ 




>o 


Therefore, our reduced Green’s function is 


(26) 


Go(x,y) = (21!) / 'Ii(k)e 


'O 


^ik-(x-y) (j){k + i/3a,x)(j){k-i/3s,y) 
F{k + il3s){\j{k + ifis) - A) 


dk. 


6. Asymptotics oe the Green’s function 

Let {ei,... ,ed) be the standard orthonormal basis in W^. Fixing u G we 

would like to show that the asymptotics (1T5|1 will hold for all {x, y) such that x — y 
belongs to a conic neighborhood containing u. Without loss of generality, suppose 
that uj ^ Cl- 

Now let IZs be the rotation in R'’* such that Tls{s) = e\ and IZg leaves the orthogonal 
complement of the subspace spanned by {s,ei} invariant. We define := 
for all j = 2,..., d. Then, (s, es,p) = (ei, Cp) = 0 and (e^,p, es,q) = (e^, e^) = 5p^q for 
p, g > 1. In other words, 

{s, 6 ^, 2 , • • •, es,d} is an orthonormal basis of R'^. 

Then around a;, we pick a compact coordinate patch so that the valued 

function e(s) = {es,i) 2 <i<d is smooth in a neighborhood of V^. 

We use the same notation for TZg and its C-linear extension to C'^. 


6.1. The asymptotics of the leading term of the Green’s fnnction. We in¬ 
troduce the function p{k,x,y) on Zi)(/co,ro) x R'^ x R'’* as follows: 


(j){k + i/3s,x)(j){k - i/3a,y) 

-■ 

where F is dehned in ([9]) and D{ko,ro) is described in Section 5.3. 

Due to Proposition 19.61 the function p is in G°°(D(/co,r) x R'’* x R'’*). For each 
{x,y), the Taylor expansion around ko of p{k) gives 

(27) p(/c, X, y) = p{ko, x, y) p'(fc, x, y){k - /cq). 
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where p' G C'°°(D(fco,I’ d) x x Note that for z & 4>{z,x) is Z'^-periodic 

in X and thus, p and p' are Z'^ x Z'^-periodic in {x,y). Since our integrals are taken 
with respect to k, it is safe to write p(fco) instead of p{ko,x,y). We often omit the 
variables x, y in p if no confusion can arise. 

Let p{k) := p{k + ko) be a cut-off function supported near 0, where p is introduced 
in Subsection 5.3. We define 


(28) 


I ■ = 

J := 


(271)--^ 

(2vr)-'^ 



^i(k-ko)-ix-y) 


^i(k-ko)-(x-y) 


dk 

Xj{k + i/3s)-X 

p{k - ko){k - ko)p'{k,x,y) 
-m-- dk 

Xj{k + ti^s) — X 


Hence, we can represent the reduced Green’s function as 


Go{x,y) = e^>^°<^-y\p{ko)I + J). 


The rest of this subsection is devoted to computing the asymptotics of the main 
integral /, which gives the leading term in asymptotic expansion of the reduced 
Green’s function Go(x, p) as \x — y\ ^ oo. 

By making the change of variables ^ = T^s{k — ko), we have 


(29) 


/ = (2i) 


—d 


P\x-y\^l_ 




-di. 


J^d {Xj o -|- TZs{ko + i/ds)) — X 

We introduce the following function defined on some neighborhood of 0 in C'^: 


Ws(z) := (A,- o n;^)i-iz + TZsiko + ifds)) - X. 


It is holomorphic near 0 (on iTZgiV)) and ITs(O) = 0. Then Wgiiz) is the analytic 
continuation to the domain TZgiV) of the denominator of the integrand in fl29l) . For 
a complex vector z = (zi,..., G C'’*, we write z = (zi, z'), where z' = {z 2 , ■ ■ ■, Zd)■ 
The following proposition provides a factorization of Wg that is crucial for com¬ 
puting the asymptotics of the integral I. 


Proposition 6.1. There exist r > 0 and e > 0 (independent of s E Vuj), such that 
Wg has the decomposition 

(30) Wg{z) = {zi-Ag{z'))Bg{z), Vz = (; 2 i,/)GH( 0 ,r)xZl'( 0 ,e)E 

Here the functions Ag, Bg are holomorphic in iA'(0,e) and 5(0, r) x iA'(0,e) respec¬ 
tively such that 24^(0) = 0 and Bg is non-vanishing on 5(0, r) x iA'(0,e). Also, these 
functions and their derivatives depend continuously on s. Moreover for z' G D'{0,e), 

(31) Ag{z') = y-Qgz' + 0{\zf), 

where 0 {\z'\^) is uniform in s when A ^ 0 and Qg is the positive definite {d — 1) x 
{d — 1) matrix 

(32) Qg = 


®See Notation 2.8 (a) in Section 2 for the definitions of i?(0, r) and D'{Q, e). 
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Proof. By Cauchy-Riemann equations for Wg and flT^ . 


OTT/" OTT/" 

(33) ^(0) = ^(0) = -7VA,(fco + 7/?.) = -VE(/3,) • s = |VE(/3,)| > 0. 


dzi 


dii 


Thus 0 is a simple zero of Wg. Due to smoothness in s of Wg and ^g, we have 


(34) 


c := min 
sev. 


^(0)> min |VE(/3,)|>0. 


Applying the Weierstrass preparation theorem (see Theorem 19.21) . we obtain the 
decomposition fl30l) on a neighborhood of 0. 

To show that this neighborhood can be chosen such that it does not depend on s, 
we have to chase down how the neighborhood is constructed in the proof of theorem 
7.5.1 in [15] (only the hrst three lines of the proof there matter) and then show that 
all steps in this construction can be done independently of s. 

In the hrst step of the construction, we need r > 0 such that Wg(zi, O') ^ 0 when 

0 < \zi\ < 2r. The mapping (s, z) ha ^ (z) = —i'VXji—iTlZ^z + ko + i/3g) ■ s is 

uZi 

jointly continuous on x IZgiV) and the value of this mapping at = 0 is greater 
or equal than c due to fl55]) and fIM]) . Therefore, ^^(z) > c/2 in some open 


dzi 


N 


neighborhood Xg x Yg of (s, 0) in x C'^. By compactness, C Xs(^k) for a 


k=l 


hnite collection of points si,..., sat on the sphere. Let Y be the intersection of all 
Ygi^ and let r > 0 such that Zl(0,2r) C Y. Note that r is independent of s. We 


claim r has the desired property. Observe that for \z\ < 2r, we have 


dWg ^ c 

dzi 2 

for any s in V^j. For a proof by contraction, suppose that there is some Zi such that 
0 < |;7i| < 2r and Wg{zi,0') = 0 = IFs(0,0') for some s. Applying Rohe’s theorem 

to the function t E [0,1] ha IFs(t2:i,0') yields ^^(fzi,0') = 0 for some t E (0,1). 

OZi 

Consequently, {tzi^Q') ^ Zl(0,2r) while \tzi\ < |zi| < 2r (contradiction!). 

For the second step of the construction, we want some 5 > 0 (independent of s) 
such that Wg{z) ^ 0 when \zi\ = r,\z'\ < 6. This can be done in a similar manner. 
Let S'(0,r) C C be the circle with radius r. Now we consider the smooth mapping 
W : {s,zi,z') HA Ws{zi.iz') where Zi E 5(0,r). Its value at each point (s,^i,0') is 
equal to IFs(2;i,0'), which is non-zero due to the choice of r in the hrst step of the 
construction. Thus, it is also non-zero in some open neighborhood Xg^zi x Yg^^i x Zg^z^ 
of (s, Zi, O') in Vuj X 5(0, r) x We select points Si,..., sm ^ and 71 ,..., '^m £ 
5(0,r) so that the union of all x Ys^^.y^,l < k < M covers the compact set 

Vi^ X 5(0,r). Next we choose 5 > 0 so that D'(0,5) is contained in the intersection 
of these Zg^^^zk- Note that 6 is independent of s and also zi. Of course IFs(^i, z') 7^ 0 


for all s and 2; E {\zi\ = r, \z'\ < 5}. According to 
in the polydisc {|^i| < r, \z'\ < 5}. 


, the decomposition 


holds 
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Also, from the proof of Theorem 7.5.1 in [15], the function Ag is dehned via the 
following formula 

(35) zi - Ag{z') = exp j ^og{zi - cu)du'j . 

The mappings (s, z') ha Ag(z') and (s, z) ha Bs(z) are jointly continuous due to fl5(l 

and fl55D . There exists 0 < e < 5 such that max|As(2;')| < r whenever \z'\ < e. We 

sev,.; 

have the identity 0301) on S(0,r) x Zl'(0,e). Now, we show that this is indeed the 
neighborhood that has the desired properties. Since \z'\ < e implies that the points 
z = { As { z '), z ') G 5(0, r) X 5'(0,e), we can evaluate (1501) at these points to obtain 

(36) Wg{Ag{z'),z')=Q, x'e5'(0,e). 


By differentiating 

(37) 


, we have 


dWg, , , dWs, , , ,x dAg, ^ r 

-^{Ag{z ),z) + -^{Ag{z ), z )-^{z ) = 0, for p = 2,..., d. 

dWg 


Observe that from the above construction, the term ^ ^ {Ag{z'),z') is always non- 

OZ\ 

zero whenever \z'\ < e. Consequently, all first-order derivatives of Ag are jointly 
continuous in (s, z). Similarly, we deduce by induction on n G that all derivatives 
of the function Ag depend continuously on s since after taking differentiation of 
the equation (13 6 p up to order n, the n-order derivative term always goes with the 
dWg 

nonzero term —— {Ag{z'), z') and the remaining terms in the sum are just lower order 

UZ\ 

derivatives. Hence the same conclusion holds for all derivatives of Bg by differentiating 


In particular, set z' = 0 in (l37p to obtain 
(38) 




dz„ 


dzi dzp 


Note that for p > 
(39) 


1 , 

dWg ^ 

OZr) 


By substituting z = 0, 

dW, 

(40) dzp 


( 0 ) 


—iVXj{—iTZg ^z + ko + ifdg) ■ TZg ^Cp. 

iV\j{ko + i/3g) ■ TZZep 

VE{Pg) ■ Cg^p = -|V5(/?s)|s • Cg^p = 0. 


(I33p . (I38l) and (I40l) imply 


dzp 


( 0 ) = 0 , 


for p 


2,...,d. 


(41) 
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Taking a partial derivative with respect to Zg {q > 1) of fl3^ at z = 0, we see that 


‘ (0) = -V(3^(*:o + iM) ■ 


dZpdZq 


(42) 


m=l 


-E 


3z 

dzm.dZ't 


'(^0 (6s,p)'m(^5, 


qjn 


m^n=l 

* Hess (yE^{^l5g^Cg^p. 

A second differentiation of flTTI) at ^ = {As{z')^z') gives 


0 = 


(43) 


(^) + -TI—(^) 


dZpdZq 


dzi 


dZpdZq 


A , d‘^Ws , .9As , d'^Ws, .dAs ,.dAs , 


dzidzq dzp 


dzpdzi dzq 


dzl dzp 


dza 


At ^ = 0, the sum in the second bracket of fl43l) is zero due to flTT|) . Thus, 


(44) 


d^As 

dZpdZq 


( 0 ) = 


dWs 

dzi 


( 0 ) 


-1 


d^Ws 

dZpdZq 


(0) {2<p,q<d). 


Together with fl3^ and 

d^A, 


(45) 


dZpdZq 


( 0 ) = 


, the above equality becomes 
1 


\VE{f3s)\2<p,,<d 


Consequently, by (HTil and (05]), the Taylor expansion of Ag at 0 implies (l3lll . 

Finally, the remainder term 0{\A\^) in the Taylor expansion (ETj), denoted by 
Rs, 3 {.z'), can be estimated as follows: 



(R Ag 

\Rs,3{^')\ ^ kT, , max 

|a|=3,0<t<l 

S.7 (*•-') 




max 
I«I=3,|i/|<|2'| 




dz^ 


(y) 


Due to the continuity of third-order derivatives of Ag on x D'(0, e), 

IRsA^') 


(46) 


lim max 

p'l^o sSVi^ 


y/|3 


< OO. 


This proves the last claim of this proposition. 


□ 


We can now let the size of the support of (d O) be small enough such that 
the decomposition (I5DD in Proposition 16.11 holds on the support of /x, i.e., supp(/i) d 
5(0, r) X D'(0,e). Therefore, from fl29p . we can represent the integral I as follows: 
(47) 


/ = (27r 


yd f ^ yy-d f f ^0 




Wg{iCj J\i'\<eJ^ 

where flg{i) = /x(C(-Bs(*0)~^- We extend jig to a smooth compactly supported 


function on 


by setting /x^ = 0 outside its support. Since all derivatives of /x^ 
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depend continuously on s, they are uniformly bounded in s. Let be the 

Fourier transform in the variable .^i of the function for each G 

i.e., 

/ + 00 
•oo 

By applying the Lebesgue Dominated Convergence Theorem, the function Ug is con¬ 
tinuous in (s, t, ^') on Vui xM"*. For such , z/s(-, i') is a Schwartz function in t on M. If 
we restrict to a compact set K containing -D'(0, e), then due to Lemma [9Tl for any 
> 0, uniformly in s and as t ^ oo. We also choose e small 

enough such that whenever |,^'| < e, the absolute value of the remainder term 0(|.^'|^) 
in fl^ is bounded from above by • QsC- Note that e is still independent of s, 
because the term is uniformly bounded by the quantity in fH 6 |) . Mean¬ 

while, each positive dehnite matrix Qg dominates the positive matrix 'Jujl(d-i)x(d-i}, 
where 7 ^^ > 0 is the smallest among all the eigenvalues of all matrices Qg (s G Vuj)- 
This implies that if 0 < |^'| < e, then 

s(i{, - A(if')) = -«(A,K')) = ■ Q,e' + o(ie7)) 

= if ■ ft?'- «(o(|??)) > > 0 . 

We thus can obtain the following integral representation for a factor in the integrand 
of / (see (HZ])): 


(48) 


- Ag{ii') 


(^ 1 , e') e M X (D'(0, e) \ {0}). 


Hence 


r-O 


/ = 


^-wAs{ii') / J{w+\x-y\)ii 


{2'kY 

1 

{2nY 


'|^'|<e J —00 
r /-h-yl 

'|5'|<e J —00 


fJ‘s{^i,^')d^idwdY 


Y-i+\^-y\)Mi^')j,^^t,Y)dtdY = h + 12 , 


where 


Ji = (27r; 




'h-2/1/2 


I 2 = (27r) 


-d 


d —00 

r rh- 2 y| 

l\^'\<eJ\x-y\/2 


e' ’'''-‘‘^‘'’>Ug{t,Y)dtdY, 


^{-t+\x-y\)Asii^') J 


Note that < x when t < \x — y\. As |a: — |/| ^ 00 , the term I 2 

decays rapidly, since for any n G N, 


( 49 ) 


hi < 


< 

rsj 



r\^-y 

'\x-y\/2 
Ax-y\ 

'\x-y\/2 


\i/g{t,Y)\dtdY 


o{\t\-^^^^^)dtdY 


o{\x-y\ "). 
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We substitute x' = \/{\x — y\ — t)^' and use the identity 


(lx -y\- t)As{ii') = --x' ■ Qsx' + O (|x'|^(|x -y\-t) 


to obtain 
(50) 


= \x-y\-(d-l)/2^^ 

^ I *'1 


-\x-y\l2 


(1 - 


exp --x' ■ Qsx' + O 


j\x'\<eyj{\x-y\-t) Y 2 V V V 

For clarity, we introduce the notation xq := \x — y\, and 


Ps t, 


/(xo,t) := 


' \x'\<e\/xo—t 


exp ( “ 2 ^^' ^ 


^Xo -1 


ps t, 


\x-y\-t 


/Xq — t 


dx'dt. 


Due to our choice of e and the dehnition of 7 ^, we have 


sup exp I — x' ■ Qsx' + O 


/Xq - t 


< exp --'yuj\x' 


Hence, for any t < xo/2, the functions 


XD'{ 0 ,e,/^t{x') ■ exp ( --x' ■ Qsx' + O 


^Xq - t 


t, 


^Xq - t 


are bounded from above by the integrable function (on 

exp {-\lu,\x'\'^\ ■ sup |z/s(t,^')l • 


Moreover, for any x' G we get 


lim XD'{o,eV^t){.x') ■ exp O 


/xo — t 


^lim^XD7o.v^)(^') ■ [t, J = P.(t,0) 

uniformly in t G (— oo,Xo/2] and in s G 

We now can apply the Lebesgue Dominated Convergence Theorem to obtain 


lim /(xo, t) = lim 


XQ^OO , 


\x'\<e^xo—t 


exp --x' ■ Qsx' + O 


, Vxo - t 

Ps(t, 0) 




= [ exp f-^x' ■ Qsx'\ Us (t, 0 ) dx' = f ey. 

J^d-i \ 2 J [det Qs A'^ 

(27r)C-i)/2 

In the third equality above, we use the change of variables u' := Qs x'. 


exp ( —] du' 


;Us{t,d). 
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Likewise, the functions X{-oo,xo/ 2 )it){l — \I{xo,t)\ are dominated (up 

to a constant factor) by the function sup|^/|<£ which is integrable on M. 

Then a second application of the Lebesgue Dominated Convergence Theorem es¬ 
tablishes 

r^o/2 

lim = lim (27r)“'^ / 


lEO^CO 


XQ^OO 

\-d 


= (2vr)- 


^ N -(d-l)/2 

1-) Hxq, t)dt 

XoJ 

foo /27r)“C-i)/2 X n 


0)dt 


hs(0) 


(27r)C-b/2(det g,)V2- 

Here we use the Fourier inversion formula for Ug in the last equality. Also, /is(0) = 

(5^(0))"^ = (^^(0))"^ = J) by dMD- Consequently, the asymptotic of the 

ozi \\/E{/3s)\ 

integral / as |a: — oo is 


I = h+l2 = 




k - si 


+ o(|l - 


(52) 


(2T)(-'-i)/2(detQ.)i/2|VB(/?.)| 


( 2 vr)C i )/2 det (-e,,p • Hess {E){^s)es,q) 2 <p,g<d 

6.2. Estimates of the integral J. 


In this part, we want to show that the expression J decays as O (|x — ■ Thus, 

taking into account ([52]), we conclude that J does not contribute to the leading term 
of the reduced Green’s function. 

In fl27|l . we set the coordinate functions of p' as (pi,... ,pd). Let us introduce the 
smooth function x, y) = pi{k + kQ, x, y)y{k) for any k G The support of 

(as a function of k for each pair {x, y)) is contained in the support of p and p^’‘\k, •, •) 
is Z'’* X ZAperiodic. We denote the components of a vector fc in M'’* as {ki,, kd). 
Observe that J is the sum of integrals Ji (1 < / < d) if we dehne 


(53) 


Ji := (27r) 


—d 


p{k-ko)-{x-y) 


y)P^^\k - ko,x,y){k - fco)/ 


'o 


Xj{k + i/3s) — A 


dk. 


Proposition 6.2. As \x — y\ —)■ oo, we have Ji = O (|a; — y\ C+i)/2^ _ 

O (lx — y\~'^/‘^) ifl>l. In particular, J = O {\x — y\~'^^‘^) ■ 


Proof. Indeed, to treat these integrals, we need to re-examine the calculation in the 
previous subsection done for the integral I. After applying the orthogonal transfor¬ 
mation IZg on each integral Ji, we rewrite them under the form of (IT7|) as 


(54) 


J, = ( 27 t)-‘ f 

J\P\<e 


^i\x-y\ii 


iii - As{if,') 


d^ide, 


where ps\^,x,y) is pd)(^^^x,y){Bs{i^)) ^ on the support of p^) and vanishes else¬ 
where. Let ,x,y) be the Fourier transform in ,^i of pi\fi,^',x,y). If the 
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parameter s is viewed as another argument of our functions here, then , x, y) 

is a Schwartz function for each quadruple It is elementary to check that 

the Fourier transform x, y) is jointly continuous on Vu) X M X ^ X X 

due to the corresponding property of j2s\^,x,y). Periodicity in {x,y) of and 
Lemma 19.11 imply the following decay: 

(55) lim |f|^ sup \ui^\t,^',x,y)\ =0, N>0. 

[s,^' ,x,y)eVu:XD'{0,e)xRd.xR<l 

In particular, 

(56) max sup \i'f\t,^',x,y)\ < oo 

(s,L|',a;,y)eV,.;XRxD'(0,e)xR'ixR'i 

and 

(57) S := mpc_ j sup ,x,y)\dt < oo. 

Jr {s,^',x,y)GVu,xD'{0,e)xR‘‘xRd' 

Recall that when 0 < |^'| < e, 3fJ(As(z^')) < 0 and thus from fl56D . 

(58) lim e^-i+\o=-y\)Asiie)^m = g. 

t^—oo 


Case 1: / = 1. 


Using 


and integration by parts, we obtain 


Ji = 


(69) 


(2ir)'' 

i 

~ {27iY 

i 

~ {2tiY 


|C'|<e J —oo 




|5'I<C J —oo 


' ^i-t+\x-y\)Min_Y^il) |/)dtde' 


dt 




^^^\\x-ylY,x,y) + 


-\x-y 


X ,x,y)dt\dY. 


Recall the notation Xq = \x — y\. The term 

[ ^i^\xo,Y,x,y)dY 




decays rapidly in xq, due to ([55]). We decompose the other term 


rXQ 


X, y)dt 

J —OO 

into two parts, where the hrst integral is taking over (xo/2,xo] and the second one 
over (— cx),Xo/2]. The hrst part decays rapidly, as I 2 in (H9|) . Now we need to prove 
that the second part decays as To do this, we use the change of variables 
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x' = \/xo — t to rewrite the remaining integral as 

fXo/2 


X, 


{d+l)/2 


' |C'|<e J —oo 


X, y)dtd(' 


(60) = 


rxo/2 


1 - 


t 

Xo 


-{d+l)/2 


' \x'\<£^X0—t 


1 / ^ / ^ / 1^' 
-X ■ QsX + O 


/|3 


X exp ( ■ Qsx' + O 


\x' 


713 


From fl57p . we derive 


(61) 


i'Xo/2 


1 - 

Xo 


^ N -(d+l)/2 


V^o -1 


sup 




X 


- t 


V^o -1, 

,x,y] dx'dt. 


Wi^\t,C,x,y)\dt < S. 


(s,^',a;,j/)sVi^xD'(0,e)xR‘^xR'^ 

On the other hand, we recall that 


3? ( --x' ■ Qsx' + O 


\x 


713 




The exponential term is majorized as follows: 


--^x' ■ Qsx' + O 


\X‘ 


713 


V^o -1 


I 1 » ^ f ^ / I X 

exp ( --X ■ QsX + O 


/|3 


\/Xo -1 


< Qx' • Qsx' + 0(e|x'H^ exp . 


Consequently, 

(62) 


'p'|<eya:o—t 


-x' ■ Qsx' + O 


i 


\x 


/|3 


V ^Xo -1 


exp ( ■ Qsx' + O 


\X' 


713 


< 

rsj 


x'l^exp ( —- 7 a;|x'|^ Idx' < oo 


\/xo -1 


7|2 \ 


dx' 


Combining (15^ through (l62|l . we deduce Ji = O(xo 


Case 2 : I > 1. 

Using fHHj) and decomposing Ji into two parts as in Case 1, we get 


r-O 


Ji = 
(63) = 


(27r)'^ 

1 

(2^ 

1 

(2^ 


/ gi(^+|x-y|)a^(/) ^7^ y)diidwdi 


|^'|<e J—oo 

r r\^-y\ 




h'|<6 J —oo 
r r\x-y\/2 

IC<6 J —OO 
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By changing the variables as before, 

r, X, y)dtde 

' J —oo 


X. 


(64) = 


|•xo/2 


1 - 


t 

Xq 


-d/2 


' \x'\<eyJxo—t 


x'l exp ■ Qsx' + ^ ^ 


\X' 


'|3 


-1 


^ I /-T 

, ^/Xo - t 


X 


, X, ?/ I dx'dt. 


In a similar manner, we obtain 


fXo/2 / , 

1 - 

Xo 


-d/2 


^x'\<e^XQi—t 


x'l exp ■ Qsx' + ^ ^ 


X 


/|3 


1 t, 


X 


yjxo - t 


.x,y 


dx'dt 


y/xo-tj ’■ 

< j |x'|exp < cx). 


This final estimate and imply J; = 0{x\ 

6.3. Estimates of the remainder term r{x,y). 


-d/2-s 


□ 


To complete the proof of Theorem 13.11 we only need to prove that the remainder 
term r(a;, y) is o{\x — y\~^f‘^^^) for any e > 0 when \x — y\ —)■ oo. The rapidly decaying 
terms I 2 and J can be ignored in the asymptotics of the remainder r, according to 
Proposition 16.21 Now let c be the coefficient in front of \x — y\~^‘^~^'>^'^ in Ji (see flS2]) ). 
Namely, c := lim where Xq = \x — y\. Hence, it is snfficient to show 

Xq^OO 

(65) lim Xo^^~^(/i — cXq = 0, 0 < £ < 1/2. 

XQ^OO 

The idea of the following lemma is to trnncate some nnnecessary (rapidly decreasing) 
parts of the main integrals we are interested in. 


Lemma 6.3. (i) For any a G (0,1) and n > 0, one has 


f‘Xo/2 


Cq J\x'\<e^y{xo-t) 


exp ( —-x' ■ Qsx' + O 


X 


/|3 


X 


t, 


X 


y/xo - t 


y/xo -t, 
dx'dt = 0(xq"'). 


(ii) Let (3 G (0,1/2), then 


'—00 JXQ<\x'\<e^(xQ-t) 


exp ( —-x' • Qsx' + O 


X 


/|3 


y/xo - t 


X 


Vs t, 


X 


y/xo -1 


dx'dt = 0{xq "■), 


for any n > 0. 
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Proof. Both parts are treated similarly to the estimates for I2. 

(i) The integral is bounded from above (up to some constant factors) by 

. ... 

x' -Qsx' ] 0 {t-^/^-^)dx'dt 


/ / ,_exp 

J x'^ J\x'\<e^{xo-t) 

xxqI2 

< / < Xr. 


(ii) We recall from fl 5 TD that on e^/xo — t), the exponential term is bounded by 
the integrable (on function 

f 1 I 

expl--7^F| 1, 

and the function 

sup |z/.(-,^')l 

is integrable over (—00,00). Hence, the integral in this lemma is estimated from 
above by: 


'—00 J XQ<\x'\<ty/^ 


exp 1 


J \2 


t, 


X 


V^o-t 


dx'dt 


< 


< 

r\j 


/ sup \us{t,C)\dt j 

' —00 sSVu,|^'|<e J 

[ exp ( — 


W \>^0 


exp ( ]dx' 




Fixing a small a, we let 

-I /»-I 

1 / 1 

/; = 


{2nr i -00 (1 - 


1 f \x' 

exp ( “2^^ ■ ^ 


/|3 


X 


t, 


X 


V^o-t 
dx'dt. 


□ 


- 1 , 

Part (i) of Lemma I 6 l 3 ] yields that Ji = Xq + 0 {xq''). Thus, we can deal with 

the truncated integral Xq instead of working directly with Ji in proving ([ 65 ]). 

So our goal is to estimate the remainder (J( — c) = 0 {xq ^^^). We represent (/( — c) 
as follows: 


( 66 ) 


exp --x’ ■ QsX' Xi-oo,x^]xD'iO,ey/^t)it,x'){l -) 


t x-(d-l)/2 


Xo' 


X exp ( O 


! |x'|^ 

V -1 


t, 


X 


_ i's{t,0)]dx'dt. 

y/xo-tj J 
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The first step is to replace the factor (l ——) in the integrand by 1. More 

xo 

explicitly, 

Lemma 6.4. IfO<a<e, we obtain 


'—oo JD'{0,e^xo—t) 


exp ( — x' ■ Qsx' + O 


\x 


/|3 




X 


yJXo - tJJ V ’ 

X ^(l-) — l^dx'dt = o{xl 

Proof. As we argued in the proof of Lemma [6.31 (ii), this integral is majorized by 


’ —OO J D' { 0 ,E\/XQ—t) 


exp( --x' • Qsx' 




X 


Vxo - t 


(l ——) — 1 ]dx'dt. 

Xo 


We break the external integral into two parts, where the first one is over (—oo, —Xq) 
and the second one is over [—Xq,Xq]. Note that the factor (l-) is less 

Xq 

than 1 if t < —Xq. Hence, the rapid decay of Og on (—oo, —Xq) (as xq is large 
enough) implies the same property for the first part. Since Ug is uniformly bounded 

on M X D'(0, e), it suffices to estimate the factor (l-) — 1 for the second 

Xq' 

part. But this is straightforward, since 


^ Xq’ 


-1 <2x“ (l-x^^) -'^-1 


= 2x“ 


1 - (1 - 


(i-x“-i)(''-b/2 
Now the claim follows from the inequality 2 q! — l<e — 1/2. 

This lemma reduces estimates the expression in fl66|) to estimating 

R^r i 


= o[x‘r^). 


□ 


exp - -x' • Qgx' ( X(-oo,x“]xD'(o,€x/S^)(^, x') exp O 


X' 


/|3 


\/Xq - t 


X Uglt, 


X 


Vxq - t 


z/s(t, 0) J dx'dt. 


Thanks to part (ii) of Lemma 16.31 it is enough to deal with the following truncation 
of i?, which can also be decomposed as the sum of and 


exp 


• Qgx^ (x(_oo,.«]xDyo,x^)(^>^')exp (o (^^== 


X Z/J t. 


X 


_, Vgit^Q'^dx'dt = R^^"^ + R^‘^\ 
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where 


= 


’ —OO j\ x'\<Xq 


1 { \x' 

exp ( — -x' ■ Qsx' + O 


/|3 


\/Xq - t 


X Z/J t, 


_ , — 0) \ dx'dt 

\/XQ-tJ J 


and 

= 


exp ( - 2 ^' •axM^(t,0)( 


X exp O 


\x' 


713 


\/xo -1 


— 1 1 dx'dt. 


Consider i?d) first. Since snp |Va;'Z^s(t, a;')| < cxo, the difference factor 

(SjpSVij xR,|a;'|<e 


I's t, 


X 


Us{t,0) 


Vxo -1 

any n > 0 large enongh, is majorized by 


is dominated by 0{\x '\/\/{xq — t)) = 0{xq For 


exp 

< «+/3-V2 

~ "''0 


-x' ■ Qsx' J dx' 


X 


/3-1/2 
0 


dt -|- 


\t\-^dt 


Now consider R^‘^\ For any /3 < 1/6, we could first estimate the below difference 
factor of the integrand 

X'(-»,.SlxD.(o, 4 )(*.^')exp (o (J=)) - 1 

< X(_oo,x«]xD'(0,x^)(^>^') (exp - l) + (1 - X(_oo,x“]xD'{0,x^)(^’^')) 

— ^(-oo,3;g‘]xD'(0,3;^)(^’ ^ )^(% ) + ^(-oo,oo)x{|a;'|>a;^} ^ ^(a:^,oo)xD'{0,x^) ^ )' 

The above inequality leads to the following estimate for : 


I A"'I < 


X, 


3/3-1/2 


|z7s(t, 0)|dt ) + / |z/s(t, 0)|dt + / exp 

J\x'\>X^ 


^0 

r*oo 


-x' ■ Qsx' ) dx' 


/ \t\-'^dt+ / r-”dr = 




Hence, R^^'> + R^‘^'> = (9(xo^^^ By choosing /? e (0, {e — a)/3), we conclude that 


R = oix, 


£-1/2n 


This proves the claim fj65ll . 


7. The full Green’s function asymptotics 


Here we show that the full Green’s function Gg^x has the same asymptotics as the 
reduced Green’s function Go as |x — 7/| —>■ cx). 
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Let US fix a point s in As we discussed in Section |5l we handle the Schwartz 
kernel Ks{x, y) of the operator Tg acting on as follows: 

r® 

= / Tg{k)d'k, 

Jo 

where JF is the Floquet transform (see Dehnition 15.ip and 

Uk) = (1 - v{k)){Ls{k) - A)-' + v{mLs{k) - \MQ,iu)))-^Qs{k). 

We would like to show that Ks{x, y) is continuous away from the diagonal and decays 
sufficiently fast as |x — |/| ^ oo uniformly with respect to s. 

Observe that the kernel of each projector Ps{k) (see Subsection 5.3) is the smooth 

function: _ 

0(fc + i/3s,x)(j){k - i/3g,y) 

F{k + i(3s) 

for each k in the support of r]. Thus, (1 — ri{k))Ps{k) is a hnite rank smoothing 
operator on T. Moreover, we also have {Ls{k) — X)Ts{k) = Ts{k){Ls{k) — A) = 
I — ri{k)Ps{k). Each Ts{k) is a parametrix (i.e., an inverse modulo a smoothing 
operator) of the elliptic operator Ls{k) — A when (s, k) G x O. This suggests to 
study parametrices of the family of elliptic operators Ls{k) — A simultaneously. 

7.1. Parameter-dependent periodic pseudodifferential operators. 

First, we briefly recall some basic dehnitions of periodic (or toroidal) pseudodif¬ 
ferential operators (i.e., TDO on the torus T). We also introduce some useful classes 
of symbols with parameters and describe some of their properties that we will use. 

There are several approaches to dehning pseudodifferential operators on the torus. 
The standard approach based on Hormander’s symbol classes (see e.g., [30]) uses 
local smooth structure on the torus T and thus ignores the group structure on T. 
An alternative approach uses Fourier series with the difference calculus and avoids 
using local coordinate charts on T (the details in Chapter 4 in [27] )@. To make a 
distinction, Ruzhansky and Turunen in [2^ refer to the symbols in the hrst approach 
as Euclidean symbols and the symbols in the latter one as toroidal symbols (see 
Section 4.5 in ra). We recall their dehnitions for only the Kohn-Nirenberg symbol 
classes, which we need here: 

Definition 7.1. Let m be a real number. 

(a) The class S'”^(T x R'^) consists of all smooth functions a{x,^) on T x R®* such 
that for any multi-indices a,(J, 

for some constant Ca^y that depends only on a, fj. Symbols in S'”^(T x R*^) are called 
Euclidean symbols of order m on T. 

(b) The class ^^(T x Z*^) consists of all functions a{x,^) on T x such that for 
each ^ G Z'^, ct(.,^) G C'°°(T) and for any multi-indices a,/3. 


different approach to periodic "ilDOs is introduced by A. Sobolev m- 
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for some constant Ca^js that depends only on a, (3. Here we recall the dehnition of 
the forward difference operator with respect to the variable ^ [27]. Let / be a 
complex-valned fnnction dehned on and 1 < j < d. Then we dehne 

A,/(0 := /(6, ■ ■ ■, 0-D 0 + 1,0+D ■ ■ ■, - /(O, 

and for any mnlti-index a, 

A“ := A“A..A“h 

Symbols in S'”^(T x Z*^) are called toroidal symbols of order m on T. 

(c) The intersection of all the classes S'”*(T x (S'”*(T x Z'^)) is denoted by 
S'“°°(T X M'^) (S'“°°(T X Z'^)), which are also called smoothing symbols. 


Dne to Theorem 4.5.3 in [27], a symbol is toroidal of order m if and only if it conld 
be extended in ^ to an Enclidean symbol of the same order m. Snch an extension is 
nniqne modnlo a smoothing symbol. Conseqnently, we will nse the notation S^(T) 
for both classes S'”^(T x and S'”^(T x Z*^). The two approaches are essentially 
eqnivalent in dehning psendodifferential operators on T whenever the symbol is in the 
class S^{T). Following [13], this motivates ns to dehne periodic psendodifferential 
operators as follows: 


Definition 7.2. Given a symbol a{x,^) G S'”^(T), we denote by Op{a) the corre¬ 
sponding periodic psendodifferential operator dehned by 

(67) {Op{a)f) (x) := ^ a{x, 


where /(^) is the Fonrier coefficient of / at The right hand side of (E7|) converges 
absolntely if, for instance, / G G°°(T). 

We also nse the notation Op(S'™(T)) for the set of all periodic psendodiherential 
operators Op{a) with a G 

Since we mnst deal with parameters s and fc, we introdnce a snitable class of 
symbols depending on parameters (s, k) G x O. 

Definition 7.3. The parameter-dependent class S'™'(T) consists of symbols (t(s, k-,x,^) 
satisfying the following conditions: 

• For each (s, k) G x O, the fnnction cr(s, fc; •, •) is a symbol in the class 

• Consider any mnlti-indices a,(3,'y. Then for each s G the fnnction 

(j(s, •; •, •) is smooth on (9 x T x , and fnrthermore, 

sup \DtDlD2iT{s,k;x,()\ < 

for some constant Cap'y > 0 that is independent of s, /c, x, and 

Tims, taking derivatives of a symbol in S'™'(T) with respect to k improves decay in 
We also denote 


mGM 


(T) := Pi 
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Definition 7.4. For each m G M U {— oo}, we denote by Op(^™(T)) the set of all 
families of periodic pseudodifferential operators {Op{a{s, k; ■, ■))}(^s,k)&d,-ixOy where 
a runs over the class S‘^{T). 

Example 7.5. 

• Suppose that |A| is small enough so that max^g§d-i |/9s| < 1. Then the family 

of symbols {(1 + {^ + k + belongs to the class ^"^(T) for any 

m G M. 

• If aa{x) G C'°°(T) and m > 0, then 

{ a„(x)(e + k + */?.)“}(.,fc) G S^{T). 

\a\<m 

• The family of elliptic operators {{Ls{k) — A)}(s^fc) is in Op{S‘^(T)). 

• If a = {a{s,k;x,^)}(s,k) e ^'(T) and b = {b{s,k;x,^)}(s,k) e ^""(T) then 
ab = {ab{s,k;x,^)}(^s,k) ^ 

• a{s,k;x,^) G S'^(T) implies D'^D^D2a{s,k;x,^) G 

The following result will be needed in the next subsection: 

Theorem 7.6. There exists a family of parametrices {As(/c)}(s,fc) in the class Op{S~‘^(T)) 
for the family of elliptic operators {{Ls{k) — A)}(s,fc)- 

The reader can refer to Section 8 for the proof of this result as well as some other 
basic properties of parameter-dependent toroidal TDOs. 


7.2. Decay of the Schwartz kernel of Tg. 


Lemma 7.7. For all k on a sufficiently small neighborhood of the support of p, A 
(< 0) is in the resolvent of the operator Ls{k)Qs{k) acting on L^(T). Furthermore, 
for such k, we have the following identity: 

(68) {{Lg{k) - A)U(Q^(fc)))-'Q,(A:) = \-^Pg{k) + {Lg{k)Qg{k) - X)-\ 

Proof. In the block-matrix form, {Ls{k)Qs{k) — A) is 

(69) 


-XPs{k) 

0 

0 

{Ls{k) - A) R(Q^(fc)) 


This gives the first claim of this lemma. The inverse of 


IS 


-X-^Pg{k) 

0 

0 

{{Lg{k) - A) R(Q^(A:)))“^ 


which proves the identity 


□ 


The identity fl68|) implies that for each (s, k), the operator 

p{k){{Ls{k) - X)\R{Q^[k)))~^Qs{k) 

is a periodic pseudodifferential operator in S~‘^{T). Thus, each of the operators Ts{k) 
is also in ^“^(T) and its symbol is smooth in (s, k) since Ps{k) and Qs{k) are smooth 
in {s,k). Actually, more information about the family of operators {Ts{k)}(^s,k) and 
their Schwartz kernels can be obtained. 
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At first, we want to introduce a class of family of operators whose kernels behave 
nicely. 

Definition 7.8. We denote by S the set consisting of families of smoothing operators 
{Us{k)}(s,k) acting on T so that the following properties hold: 

• For any mi, m 2 G M, the operator Us{k) is smooth in A; as a H^^{T))- 

valued functiorQ. 

• The following uniform condition holds for any multi-index a: 

s,k 

We remark that if the family of smoothing operators {Us{k)}(^s^k) is in 
then this family also belongs to S. 

In order to obtain information on Schwartz kernels of a family of operators in S, we 
need to use the following standard lemma on Schwartz kernels of integral operators 
acting on T. 

Lemma 7.9. Let A be a hounded operator in L^(T). Suppose that the range of A is 
contained in where m > d/2 and in addition, 

for all f 

Then A is an integral operator whose kernel KA{x,y) is bounded and uniformly 
continuous on T x T and the following estimate holds: 

(70) lKA(x,g)l < joC, 

where 70 is a constant depending only on d and m. 

The fact can be found in Lemma 2.2 in [2]. 

Now we can state a useful property of Schwartz kernels of a family of operators in 

5 . 

Corollary 7.10. If {Us(k)}(^s,k) is a family of smoothing operators in S, then the 
Schwartz kernel Ku^{k,x,y) of the operator Us(k) satisfies 

sup \D'/Ku,{k,x,y)\ < 00 , 

s^k,x,y 

for any multi-index a. 

Proof. We pick any m > d/2. Then by Dehnition 17.81 we have 

sup \\D//Us(k)f\\H”^(T) < Ca\\f\\H-^(T)- 

s,k 

Applying Lemma 17^ the estimates flTOj) hold for kernels D'^Ku„(k,x,y) of the oper¬ 
ators D'^Us(k) uniformly in (s, k). □ 

We now go back to the family of operators Tgik). 

^We remind the reader that B{E,F) denotes the space of all bounded linear operators from the 
Banach space E to E. 





SPECTRAL EDGE INTERIOR CASE 


33 


Proposition 7.11. There is a family of periodic pseudodifferential operators {Bs{k)}(s,k) 
in Op{S~'^{T)) such that the family of operators {Ts{k) — Bs{k)}{^s,k) belongs to S. 

Proof. Due to Theorem 17.61 there is a family of operators in Op{S~‘^(T)) 

and a family of operators in such that 

{Ls{k) - X)As{k) = I - Rs{k). 

Since Ts{k){Ls{k) — A) = / — ri{k)Ps{k), we deduce that 

(71) Tsik) = Mk) - r]{k)Psik)Asik) + Tsik)Rsik). 

One can check that the symbols of operators ri{k)Ps{k)As{k) are in the class 

due to the composition formula. Thus if we let Bs{k) := As{k)—T]{k)Ps{k)As{k), 
then {Bs{k)}(^,^k) e Op{S-‘^{T)). 

Now it remains to show that the family of smoothing operators {Ts{k)Rs{k)} 
is in S. Let us hx any two real numbers mi, m 2 and a multi-index a. Notice that 
{Ls{k) — A) is analytic in /c as a i7™'^“^(T))-valued function and also, 

snp \\D^{Ls{k) — < 00 . 

s^k 

Due to Lemma 17771 

Tg{k) = (1 - v{k)){Lg{k) - A)-i + rj{k)\-^Pg{k) + p{k){Lg{k)Qg{k) - \)-\ 

Thus, Ts{k) is smooth in fc as a iL™'2(^']r))_valued function and moreover, 

(72) snp\\D^Ts{k)\\B{H^ 2 -^{T),H^ 2 {T)) < 00 . 

s^k 

Since {Rs{k)} is in Op(^“°°(T)), Rs{k) is smooth in A; as a B(T), (T))- 

valued function and furthermore, 

(73) sup ||D^i?s(A:)||5(l^mi(T) l^m 2 - 2 (T)) < 00 . 

s,k 

By fl72|) . fl73l) and Leibnitz’s rule, we deduce that Ts{k)Rs{k) is smooth in fc as a 
B(/7"*i(T), i7”^2(T))-valued function and the corresponding uniform estimate also 
holds. We conclude that the family {Ts{k) — Bs{k)}(^s,k) is in S. □ 

We need the following important estimate of Schwartz kernels of operators Ts{k): 

Corollary 7.12. Let Ks{k,x,y) be the Schwartz kernel of the operator Ts{k). Let 
N > d — 2. If a is a multi-index such that |a| = N, then each D'^Ks{k,x,y) is a 
continuous function on T x T and the following estimate also holds uniformly with 
respect to {x, y): 

sup \DfKs{k.iX^y)\ < 00 . 

(s,fc)esrf-ixC) 

Proof. Due to Proposition 17.111 the operator Ts{k) is a sum of operators Bg{k) and 
Ug{k) such that {i?s(A;)}(s,fc) € Op(^“^(T)) and {Ug{k)}(^g^k) G S. In particular, 

Kg{k,x,y) = KB,{k,x,y) + Ku,{k,x,y). 
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Recall that in the distributional sense, the Schwartz kernel KB^{k,x,y) of the 
periodic pseudodifferential operator Bs{k) is given by 

where (t(s, k;x,^) is the symbol of the operator Bs{k). 

Since is in 5“^(T), 

< (1 + 

Since — (2 + iV) < —d, the sum 

converges absolutely and moreover, 

sup \DtKB.(k,x,y)\ < V(1 + lel)-'"" < ex., 

{s,k,x,y)£S‘^ ixOxTxT 

Combining this with Corollary 17.101 we complete the proof. □ 

Notation 7.13. Let be a function on and 7 be a vector in M^, then T^'ip is the 
7 -shifted version of "0. Namely, it is defined as follows: 

+ 7)- 

We denote by V the subset of consisting of all functions -0 such that its 

support is connected, and if 7 is a non-zero vector in Z'^, then the support of r^yip 
does not intersect with the support of 0 . 

Definition 7.14. Since is the universal covering space of T, we can consider the 
covering map 

TT : ^ R'^/Z'^ = T. 

In particular, ■n[x + 7 ) = vr(a;) for any a: G R'^ and 7 G Z*^. 

A standard fundamental domain (with respect to the covering map vr) is of 
the form [0, l]*^ + 7 for some vector 7 in R'^. Thus, a standard fundamental domain 
is a fundamental domain of R'^ with respect to the lattice Z'^. 

Using Definition 15.11 of the Floquet transform we can obtain the following 
formula: 

Lemma 7.15. Let 0 and 9 he any two smooth functions in V. Then the Schwartz 
kernel of the operator (fTgO satisfies the following identity for any (x, y) G 

R'^ X R*^.- 

Ks,^,eix,y) = [ e''^'^^~y^(t){x)Ks{k,T:{x),'K{y))e{y)dk. 

Jo 

Proof. Since both 0, d G V, there are standard fundamental domains lU^ and Wy C 
R'^ so that 


supp(0) C supp(6') C We. 
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Then, it suffices to show that {(^TgOf, g) equals 

1 r f 


(2i) 


'W^ JWe Jo 


y){x)Ks{k, 'k{x ), 7i{y)){9f){y)dkdydx, 


for any /, g in 
We observe that 


mef,g) = {J^(j)T,ef,Tg) 
1 


{27rY 

1 


(^£ T,{k)dk^ H0f ),^9 


'o 


T,{k)dk)j^i9f),j^{(j)g) 


\ 




{2tiY 

Since 9 eV, for any y in Wy, we have 

:F{9f){kMy)) = {0f){y)e-^^-y. 

Similarly, 

J^((i>g){k,Ti{x)) = {(i>g){x)e~'^^'^, Vx e W^. 

We also have 


TYk)dk]{nomk,n{x)) = TYk){HOf){kr)){Ax)). 


'o 


/ 


Consequently, 


'o 


TYk)dk]J^{9f),J^{cl>g)'j 



Ts{k){T{9f){k, ■)){TT{x)){(j)g){x)e~'^^'^dxdk 

'o Jw^ 

III Ks{k,TT{x),TT{y))J'{9f){k,TT{y)){(j)g){x)e'''''''dydxdk 
'o Jw^ JWe 

f e^H^~y)Ks{k, Ti{x) , TT{y)){9f)(y){(jyg){x)dydxdk. 

>0 JW 4 , JWe 

Using Fubini’s theorem to rewrite the above integral, we have the desired identity. □ 



Proposition 7.16. Consider any two smooth compactly supported functions 0 and 9 
on Mf such that their supports are disjoint. Then the kernel Kg^^^x, y) is continuous 
on X and moreover, it satisfies the following decay: 

sup \Ks^^Yx,y)\ < CN\(t){,x)9{y) \ ■ |a; - y\~^, 

S 

for any N > d — 2. Here, the constant is independent of fi and 9. 


Proof. By using partitions of unity, any smooth compactly supported function can 
be written as a hnite sum of smooth functions in the set V. Thus, we can assume 
without loss of generality that both 0 and 9 belong to V. 

First, observe that for any {k,n) E O x ifi, 

Tfik + 2nn) = MffiTs{k)Mn, 
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where Mn is the multiplication operator on L^(T) by the exponential function 
Hence, 

V^K,(k + 27 m, 7 r(a:), 7 r(y)) = 7 r(a:), 

for any multi-index a. Since for any x E M'^, we obtain 

(74^ ^i(k+2nnHx-y)yaj^^^^ ^ 27m, 7r(a:), 7r(y)) = 7r(x),7r(y)). 

Applying Lemma lY.151 we then use integration by parts (all boundary terms vanish 
when applying integration by parts due to (|74|l ) to derive that for any |a| = JV, 

{2TTf{i{x-y)YKs,<l,fi{x,y) = 4>{x)e{y) [ e^^'^^~^^VkKsik,Tr{x),TT{y))dk. 

Jo 

Suppose N > d — 2. Then by applying Corollary 17.121 the above integral is 
absolutely convergent and it is also uniformly bounded in (s, x, y). Consequently, the 
kernel Ks^,pfi{x, y) is continuous. Furthermore, 

sup I a:,, 0 , 0 ( 0 ;, 7 /) I < \(f){x)e{y)\ - min |(o:-7/)“|"^ < |0(x)0(7/)| • |x-7/|"^. 

s \a\=N 

□ 

We now have enough tools to approach our goal: 

Theorem 7.17. The Schwartz kernel Ks{x, y) is continuous away from the diagonal 
and furthermore, as \x — y\ ^ 00 , we have 

sup\Ks{x,y)\ = 0{\x-y\~^), ViV > 0. 

S 

Proof. Let us £x a point (s, x) in x W^. Now we consider a point y = x + st, 
where f is a real number. When |f| > 0, we can choose two cut-off functions (j) and 
6 such that 0 and 6 equal 1 on some neighborhoods of x and y, respectively, and 
also, the supports of these two functions are disjoint. Then, Proposition 17.161 implies 
that the kernel Ks{x,y) is continuous at {x,y) since it coincides with on a 

neighborhood of {x,y). This yields the first claim. Again, by Proposition 17.161 we 
obtain 

sup \Ks{x,y) \ = sup |A'^,0,0 (x,7 /)| < Cn\x - y\~’^, 

S S 

which proves the last claim. □ 

8. Some results on parameter-dependent toroidal TDOs 

The aim in this section is to provide some results needed to complete the proof of 
Theorem 17.61 We adopt the approach of [13] to periodic elliptic differential operators. 

The next two theorems are straightforward modihcations of the proofs for non¬ 
parameter toroidal TDOs: 

Theorem 8.1. (The asymptotic summation theorem) Given families of sym¬ 
bols bi E S”*“^(T), where each family bi = {6/(s, fc)}(s,fc) for I = 0,1,..., there exists 
a family of symbols b in S^(T) such that 

(75) {b{s, k)-J2 Hs, k)}is,k) e ^"^-'(T). 

i<l 
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We will write b rsj E bi ifb satisfies m 

I 

Proof. Step 1. Let n = m + e for some e > 0. Then 
\b,(s,k;x,()\< c,(i + lei)”-' = 

Thus, there is a sequence {rii}i>i such that rji +cxo and 

for |,^| > Tji. Let p e C°°(M) satisfy that 0 < p < 1, pit) = 0 whenever |t| < 1 and 
pit) = 1 whenever |t| > 2. We dehne: 


b{s,k-x,f) = ^p{^^^ 


bi{s,k-,x,f). 


Since only a finite number of summands are non-zero on any compact subset of TxM'^, 
b{s, •; •, •) G C°°{0 X T X M'^). Moreover, b{s, k) — J2r<i br{s, k) is equal to: 


X] ( P ( — ) - 1) br{s,k) + bi{s,k) + ( — ) br{s,k). 


Pr 


lei 


Pr 


r<l \ \ 

The first summand is compactly supported while the second summand is in 
Now let e < 1. Then, the third summand is bounded from above by 




m—l 


r>l 


Consequently, 


sup 


6 (s, k) — br{s, k) 


r<l 


<^^(l+lel) 


m—l 


Step 2. For |a| -|- |/3| -|- I 7 I < N, one can choose pi such that 


sup \DfD^Dfbi{s,k;x,f)\<^{l + \f\) 

Sg§d-1 / 


n-l-\a\-\0\ 


for p; < 1^1. The same argument as in Step 1 implies that 


(76) sup \DlDlD-<Ms,k) ~y^K(s,m < C„(l + 

ses<i-i 

Step 3. The sequence of pfs in Step 2 depends on N. We denote this sequence by 
Pi^jq to indicate this dependence on N. By induction, we can assume that for all I, 
Pi,N < Pi,N+i- Applying the Cantor diagonal process to this family of sequences, i.e., 
let Pi = pi^i then b has the property 0761) for every N. □ 


Theorem 8.2. (The composition formula) Let a = {a(s, P)} be a family of 
symbols in ^*(T) and P{x,D) = J2\a\<m^a{^)D^ be a differential operators of order 
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m > 0 with smooth coefficients affix). Then the family of periodic pseudodifferential 
operators {P{xffi + k + i/3s)0p{a{s, k))}(^s,k) ^ Op{S’‘~^"^(T)). Indeed, we have: 

Pffk)Op{a{s, k)) = Op{P o a)(s, k)), 


where 

(77) 

and 


{Poa){s,k]xffi) = —DlPffk){xffi)D^a{s,k]xffi) 


|/3|<m 


Ps{k){xffi) = P{xffi + k + i/3s). 


Proof. The composition formula fl77|) is obtained for each (s, k) is standard in pseu¬ 
dodifferential operator theory (see e.g., [T3l[23[30]). We only need to check that the 
family of symbols {(P o a)(s, k;x,f)}(^s,k) is in S'^+™(T). But this fact follows easily 
from (1771) and Leibnitz’s formula. □ 


We now hnish the proof of Theorem 17.61 


Theorem 8.3. (The inversion formula) There exists a family of symbols a = 
{a(s, in S'“^(T) and a family of symbols r = {r(s, in S~°°(T) such 

that 


{Lsik) - A)Op(a(s, k)) = I - Op{r{s, k)). 


Proof. Let 


and 


Lo{s,k;xffi) := Y affix){f + k + 

\a\=2 

||o||oo := 'Y II®«(‘)IU°°(t)5 

l«l=2 


M := max (\k\‘^ + 9 ^||a|loo|/5sP + 6^ , 

(s,fc)6Srf-ixC) ^ ^ 


where 9 is the ellipticity constant in ([2]). Whenever |^| > 


\Lo{s,k-,xffi)\ > 3f?(Lo(s, A:;a;,0) > + k\^ - Y 

l«l =2 

>9(W-\k\^^ - ||a|UI/5.r > 1. 

Let p e C°°(R) be a function satisfying p{t) = 0 when \t\ < and p{t) 

when |f| > 2M^/^. We dehne the function 


(78) 


ao{s,k){xffi) = p(l^l) 


1 

Lffis,k;xffi)' 


1 


Then oq := {ao(s, A:)}(s,a;) is well-dehned and belongs to S ^(T). The next lemma is 
the hnal piece we need to complete the proof of the theorem. 
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Lemma 8.4. (i) If b = {b{s, k)}(^s,k) G S’'{T) then b — {L — X) o [aob) G 

(ii) There exists a sequence of families of symbols ai = {a;(s, k)}(^s^k) in I = 

0,1,... and a sequence of families of symbols ri = {ri(s, A:)}(s,fc) in S~\T), Z = 0,1,... 
such that qq is the family of symbols in m. ro(s, k) = 1 for every (s, k) and for all 

I, 

{L - \) o ai = ri - ri+i. 

Proof (i) Let p{s,k) = {L{s,k) - A)(x,0 - Lo{s,k-,x,^) so that p = {p{s, k)}(^s,k) e 
S'^(T) and hence, p o (aob) is in due to Theorem 18.21 Moreover, b — Loaob = 

(1 — p(|^|))& is a family of symbols whose ^-supports are compact and thus it is in 
We can now derive again from the composition formula 0771) when P := Lq 

that 

(L - A) o (aob) = Lqo (aob) +po (aob) = Loaob H-= 6 + ..., 

where the dots are the terms in S'‘~^(T). 

(ii) Recursively, let ai = aoVi and r;+i = ri — (L — X) o ai. By part (i), r;_|_i e 

Now let a be the asymptotic sum of the families of symbols ai, i.e., a 
Then 

(L - A) o a ~ - X) o ai = '^ri - n+i = ro = 1, 

i i 

which implies that 1 — (L — A) o a ~ 0. In other words, this means that r : = 
1 — (L — A) o a G S~°°(T). Hence, there exists a family of symbols a in S~‘^(T) and a 
family of symbols r in S~°°(T) satisfying (L — A) o a = 1 — r. Finally, an application 
of Theorem 18.21 completes the proof of Theorem 17.61 □ 


9. Some auxiliary statements 

9.1. A lemma on the principle of non-stationary phase. 

Lemma 9.1. Let M be a compact manifold (with or without boundary) and a : 
R X M —)• C &e a smooth function with compact support. Then for any N > 0, there 
exists a constant Cn > 0 so that the following estimate holds for any non-zero t G R.' 


(79) 


sup 

xGM 


e"^^a(y,x)dy 


< 


-N 


Here depends only on N, the diameter R of the y-support of a and sup {d); a 




Proof. Let t 0. Applying integration by parts repeatedly (A-times), it follows that 


e**^a(|/, x)dy 


= \t\ 


-N 


P^yvy(y,x)dy 


< Rsup \dy a\ ■ \t\ 


x,y 


□ 
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9.2. The Weierstrass preparation theorem. 

Theorem 9.2. Let f{t,z) be an analytic function of{t,z) G in a neighborhood 
o/(0,0) such that (0,0) is a simple zero of f, i.e.: 

/( 0 , 0 )= 0 , ^( 0 , 0 ) ^ 0 . 

Then there is a unique factorization 

f{t,z) = {t- A{z))B{t,z), 

where A, B are analytic in a neighborhood of 0 and (0, 0) respectively. Moreover, 
5(0,0) ^ 0 andA{0) = 0. 


The proof of a more general version of this theorem could be found in Theorem 

7.5.1 in [I5]. 

9.3. Proofs of Proposition 4.1 and Lemma 5.3. 

Remark 9.3. Consider a domain B of and let / : P —)■ C be a holomorphic 
function. For 2 ; G C'^, write z = x + iy where x,y E Now we £x a vector (3 in 
and denote = ifD — ifd) fl M'^. If this intersection is non-empty, we may consider 
the restriction k -E- f{k + i/3) as a real analytic function dehned on a subdomain 
Bjs of M'^. Thanks to Cauchy-Riemann equations of /, we do not need to make any 
distinction between derivatives of / with respect to x (when / is viewed as a real 
analytic one) or z (when / is considered as a complex analytic one) at every point in 
Vp since 

df df df 

— (fc i/3) = -TT—ik + i/3) = —i—{k i/3), 1 < I < d. 

0 x 1 ozi oyi 

For higher order derivatives, we use induction and the above identity to obtain 

crj(k + i/3) = &;f(k + iD) = (-i)\-\dkp(k + !/3), 

for any multi-index a. We use these facts implicitly for the function Aj. When dealing 
with the analytic function / = Aj in this part, denote d"Aj to indicate either its x or 
^-derivatives. 

We also want to mention this simple relation between derivatives of \j and E\ 

a“5(/3) = a“A,(fco + 1 / 3 ) = ^i“ia“A,(fco + */?). 

Proof of Proposition f.l. We recall from Section 2 that V is an open neighborhood 
of ko in such that the properties (P1)-(P6) are satisfied. Note that V depends 
only on the local structure at fco of fhe dispersion branch Xj of L. Denote Og = 
{k + itf3s : k e 0,t E [0,1]} for each s G For C > 0 (which we define it later), 

set Mg^c = Og n {z E : |3?(2;) — /co| < C} and W,c = \ Ms,c- For C and |A| 

small enough, we can suppose Mg^c V since /3g is small too. We also assume that 

|A| < eo@ 

For any point z = k + itf3g E Mg^, we want to show if A G a{L{z)), it forces 
z = ko i(3s. By (P3), this is the same as showing the equation Xj{z) = A has 
no solution 2 ; in Mg^c except for the trivial solution z = ko + i/3g. Suppose for 


’Recall the definition of eg from (P3). 
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Figure 1. An illustration of the regions Mg^c and Ng^c when = 0. 


contradiction Xj{k + it(3s) = X = X{(3s) for some t G [0,1] and A; in {fc G | 0 < 
\k — ko\ < C}. By Taylor expanding around fco + At/ds, there is some 7 G (0,1) such 
that 


(80) 

A - ^ 


{k - ko) ■ VA,(fco + itPs) + + It^ 


a\ 


|a|=3 

{k - koY 


+ -{k -ko) ■ Hess {Xj){ko + it/3g){k - ko) + - -^d°'Xj{'y{k - ko) + ko + it^s) 

\ z ^' al 


\a\=A 


If |a| is odd, then by Remark 19.31 and the fact that E is real, we have 

a"Aj(fco + itl3g) = -^^d°'E{tl3g) G *. 

Taking the real part of the equation flSU]) to get 

E{l3g) - E{tl3g) = -hk - ko) ■ Hess(F)(t/3,)(fc - fco) + V ~ x 
/ ^^ al 

\a\=A 

X ^{d°‘Xj{'y{k - ko) + ko + it/3g)). 


The left-hand side is bounded above by (1 — t)A < 0 because of the concavity of E 
{E(t(5g) > tE{l3g) = tX). On the other hand, by (P5), 

--(A) — ko) ■ Hess {E){t/3s){k — ko) > - \k — AjqP mincr(Hess (Aj)(0)) 


E ^{d^Xj{j{k - ko) + ko + it(3g)) 


al 


|a|=4 

We simply choose < 


2 ^ min (T(Hess (Aj)(0)) 


C{d) max^g^ |^|^ja“Aj(z)| 


< C{d)\k — ko\‘^ max \d°‘Xj{z)\. 

zeV,\a\=4 

to get a contradiction if A: 7 ^ fco- 
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For the remaining part, we just need to treat points k + it 13s in Ng^c- We have 
A G p{L{k)),Wk e m3. The idea is to adapt the upper-semicontinuity of the spectrum 
of an analytic family of type A on following na. For any k ^ O and z G C®*, the 
composed operators {L{k + z) — L{k)){L{k) — A)“^ are closed and dehned on L^(T) 
and by closed graph theorem, these are bounded operators. Clearly, 

L{k + z)-\ = {l + {L{k + z)- L{k)){L{k) - X)-^){L{k) - A). 

Thus, A is in the resolvent of L{k + z) if the operator 1 + {L{k + z) — L{k)){L{k) — X)~^ 
is invertible. Hence, it is enough to show that there is some positive constant r such 
that for any k E O and \z\ < r, 

(81) \\{L{k + z)- L{k)){L{k) - A)-i„p < 1/2, \k - ko\ > C, 

where the operator norm on L^(T) is denoted by || • ||op. Indeed, if |A| is small enough 

so that we have max |/3s| < r and then fISTD implies that A G p{L{k + it(3a) for any 
seS'^-i 

fG [ 0 , 1 ]. 

Finally, we will use some energy estimates of linear elliptic equations and spectral 
theory to obtain flST]) . Observe that, 

L{k + z) — L{k) = • A{x){D + k) + {D + k) ■ A{x)z + z ■ A{x)z. 

For V G H^(T) and \z\ < 1, there is some constant Ci > 0 (independent of z) such 
that 


(82) \\{z ■ A{x){D + k) + {D + k) ■ A{x)z + z ■ A{x)z)v\\l 2 (t) < Ci\z\ ■ ||u||hi(t)- 

Set V := {L{k) — X)~^u for u G L^(T). Ellipticity of L(k) yields v G H‘^{T) and in 
particular, we obtain flH^ for such v. Testing the equation (T(/c) — A)u = u with the 
function u, we derive the standard energy estimate 


(83) ||Zlu||i2(']r) < C2(||u||2,2(']r) -1- ||u||i2(']r)). 

Note that both Ci and C 2 in flH^ and fl8^ are independent of k and A since we take 
k in the bounded set O and consider |A| to be small enough. 

Suppose that |A| is less than one-half of the length of the gap between the dispersion 
branches Xj and Aj_i. Due to functional calculus of the self-adjoint operator L(/c), 
we get 


||(L(fc)-A) ^\\op = dist{X,a{L{k))) ^ = min{(Aj(fc) - A), (A - Aj_i(fc))} h 


Now let = —|maxAj_i(fc) > 0 and 82 
(^2 > 0. Moreover, 


min Xj(k). Then due to A3, 
fcec>,|fc-fcol>c 


A — Aj_i(/c) > A 


max Xi-i(k) > (5i, 
fceo 


and 


Xi(k) — A > min AHfc) — X > 62 . 

^ - k&0,\k-ko\>C ^ 


Hence, 

( 84 ) 
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In other words, ||T||L 2 (']r) < ( 5 ||M||i 2 (T). Applying this fact together with fl 8 ^ and 
dSS]), we have 

||(L(fc + z) — L{k)){L{k) — A) ^u\\l 2 (j^ < IzIC'iUnH/fpT) 

< |2|CiC2(||t||l2(t) + ||'?^IU2(t)) 

< \z\CiC2{)~ + 5)\\u\\i2[jy 

Now flSip is a consequence of the above estimate if we let 

□ 

Proof of Lemma 5.3. From Lemma I2.6[ the complex Bloch variety E := Bl of the 
operator L is an analytic subset of codimension one in By [I9l|33], there exist 

an entire scalar function h{k,fi) and an entire operator-valued function I{k,fi) on 
such that 

( 1 ) h vanishes only on S and has simple zeros at all points of S. 

(2) In \ S, {L{k) - /i)-i = h{k, fi)-^I{k, /i). 

In particular, {Lt^s{k) — A)“^ = h{k -f itjSs, X)~^I{k + it/3s,X) for fc G and 
t E [0,1) by Proposition |4.11 Due to Assumption A and (P2), ifi„ -|- it/3g E Id, 
the fc-variable function h{k,X)~^ is equal (up to a non-vanishing analytic factor) to 
{Xj{k + it/3s) — A)“^ on an open disc D{ko, 2e) C 1/ in C'^ for some e > 0. Hence, we 
can write the sesquilinear form for such values of k as 


{Rt,s,xf, p) — Ri + R2) 


where 


Ri 


JonD{ko,e) 


\j{k + it/Ss) — A 


and 

R2 = [ (L{k + it^g) - xy^fik), yk)) dk. 

Jo\D{ko,e) ^ 2 

Here M(fc, A) is a L^(T)-valued analytic function on D{kQ, e) when |A| is small. Since 
/ and ip have compact supports, their Floquet transforms f{k), (p{k) are analytic with 
respect to k. To prove the equality fl2^ . we apply the Lebesgue Dominated Conver¬ 
gence Theorem. For Ri, it suffices to show that the denominator in the integrand 
when t —)• 1“ is integrable over D{ko,e) for d>2. Indeed, 


\Xj{k + i(3s) - A| > 5 


zVE{y) • {k - ko) 


ko) ■ Hess {E){f3s){k - ko) , 


for some 5 > 0 if e is chosen small enough so that in the Taylor expansion of Xj at 
ko + i(3s, the remainder term 0(|fc —/cqI^) is dominated by the quadratic term \k — ko\^. 
Furthermore, 


zVE{y)-\k-ko 


^{k- ko) ■ Hess {E){/3s){k 


ko) 


> C{\k - ko\^ + \k - kol^), 
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for some constant C > 0 (independent of k). Now let v := {k — k^) and so the right 
hand side of the above estimate is just |np + (up to a constant factor). One can 
apply Holder’s inequality to obtain 

> C\vif’/‘^\v'\, 

where v = (ni,n') G M x Thus, we deduce 

(85) \Xj{k + i^s) - Al"^ < C\vi\-^/^\v'\-^/^ 

Since the function \x\~^ is integrable near 0 in M'’* if and only if n > d, and 

are integrable near 0 in and M respectively. Therefore, the function in 
the right hand side of fl55]l is integrable near 0 . 

The integrability of i ?2 as t —)■ 1“ follows from the estimation fISTD in the proof of 
Proposition 14.11 Indeed, 

(86) 

\\{L{k + ztp,) - A)-lop = 11(1 - {L{k + zt/3s) - L{k)){X - L{k))-Y\X - 

<_ mk)-X)-Xp _ 

- 1 - \\{L{k + itX - L{k)){X - L{k))-Xp' 

By decreasing |A|, if necessary, and repeating the arguments when showing flHT]) and 
flHTll we derive: 

(87) 1 - \\{L{k + tt/3s) - L{k)){X - L{k)XX > 1/2, Vfc G O \ D{ko,e) 

and 


( 88 ) 


Thanks to 


\-ii 


sup ||(L(/c)-A)' nop 

keO\D{ko,6) 


< oo. 


, (157)) . (1551) . Cauchy-Schwarz inequality and Lemma 15^ we have: 


sup 

te[o,i] 


and 


L{k + tt/3s) - XXf{k), ip{k) < 2|| {L{k) - A) 


\-ii 


I 


||/(fc)||L2(T)||<7’(fc)||L2(T) 


JO\D(ko,E) 

This completes the proof of our lemma. 


<\\fmM\m\XT) ,ykeo\D{ko,e) 

\\fX\\L^{T)\Mk)\\L^T)dk < ||/||L2(Rd)||(p||i2(Rd) < CX). 


□ 


9.4. Regularity of eigenfunctions (j){z,x). 

In this subsection, we study the regularity of the eigenfunctions (j){z, x) of the 
operator L{z) with corresponding eigenvalue Xj{z) (see (P4)). It is known that for 
each z G V, the eigenfunction (j){z, x) is smooth in x. We will claim that these 
eigenfunctions are smooth in {z,x) when 2 ; is near to fco- The idea is that initially, 
(l){z, ■) is an analytic section of the Hilbert bundle V x H‘^{T) and then by ellipticity, 
it is also an analytic section of the bundle V x H'^{T) for any m > 0 (for statements 
related to Fredholm morphisms between analytic Banach bundles, see e.g., [35]) and 
hence smoothness will follow. 
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For the sake of completeness, we provide the proof of the above claim by applying 
standard bootstrap arguments in the theory of elliptic differential equations. 

Lemma 9.4. The function df(j){z,x) is jointly continuous on 1/ x for any multi¬ 
index a. 


Proof. By periodicity, it suffices to restrict x to T. Let K := V. Due to (P4), the 
function z i—)■ •) is a iL^(T)-valued analytic on some neighborhood of K. Thus, 

sup \\(j){z, •)||ri2(T) < oo. 

z&K 

Then, we can apply bootstrap arguments for the equation 

L(z)0(z, •) = \j{z)(t){z,-) 

to see that Mm := sup^gic ||0(^, is finite for any nonnegative integer m. 

Now we consider 2 ; and z' in K. Let (fz^zfx) := (f){z,x) — (j){z',x). Then, (j)z,z' is a 
(classical) solution of the equation 

L{z)(j)z,z' = /z.P, 

where := {\j{z)(t){z, ■) - \j{z')(t){zj •)) + {L{z') - L{z))(t){zj •)• 

By induction, we will show that for any m > 0, 

( 89 ) \\(t^z,z'\\H-^{T) ^\z — z'\. 

The case m = 0 is clear because (P4) implies that 2 ; 1 —)■ ||0(2 ;, OIUfT i® Lipschitz 
continuous. 

Next, we assume that the estimate fl8^ holds for m. As in fl8^ . 

(90) ||(L( 2 ;) - L{z'))(l){zj Olk-m < k - z'\ ■ 110 (2;', < Mm+i\z - z'\. 

Using triangle inequalities, the estimates fl89D . fl90|) and analyticity of Aj, we get 

(91) 

\\fz,z'\\H^{T) ^ \\^j{z)4>{z, ■) - Xj{z')4>{z', •)||h"‘(T) + \\{L{z) - L{z'))4>{z', •)||h"‘(T) 

^ ■ || 02 , 2 '||ir'"(T) + Mm\Xj{z) — Aj( 2 ;')| + Mm+i\z — z'\ 

<\z-z'\. 

Notice that for any m > 0, the following standard energy estimate holds (see e.g., 

P[niE3|): 

(92) ||02,2'||h’"+2(t) ^ \\fz,z'\\H^{T:) + 

Combining fl^ and fl9^ . we deduce that \\f>z,z'\\H^+'^{r) \z — z'\. Hence, flS^ holds 
for m + 2. This finishes our induction. 

Applying the Sobolev embedding theorem, ||0 ^,^'||c™(t) '^\z — z'\ for any m > 0. 
In other words, 0 G C{K.,C'^(T)) for any m. Since C{K x T) = C{K,C(T)), this 
completes the proof. □ 


Notation 9.5. Consider a z-parameter family of linear partial differential operators 
{L( 2 :)} where 2 : G Suppose L{x,f,z) is the symbol of L{z). Whenever it makes 

dLiz) dL 

sense, the differential operator ——- is the one whose symbol is -^{x,f,z) for any 

OZi OZi 

/ G {1, 2, ..., d}. 
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Proposition 9.6. Assume D is an open disc centered at ko in such that D ± 
ifds (E V for any s G Then all eigenfunctions (j){k ± i/3s,x) are smooth on a 

neighborhood of D x Furthermore, all derivatives of (j){k ±il3s,x) are bounded on 

D X R'^ uniformly in s, i.e., for any multi-indices a, 13: 

sup \d'^d^(l){k ±if3s,x)\ < oo. 

{s,k,x)£S‘^~^ X DxR‘^ 


Proof. Pick any open disc D' in R'^ so that D ± ijdg C -D' ± ijdg C V. We will prove 
that all eigenfunctions are smooth on the domain D' x R*^. Also, it is enough to 
consider the function (f(k + i/d^) since the other one is treated similarly. 

d(j) 

First, we show that Tvrik + il3s,x) is continuous for any 1 < / < d. By Lemma 
oki 

I9.4[ the function (fc, x) i—)■ (fik-\-i(3s, x) is continuous on D' x T. We consider any two 
complex-valued test functions tf G C^(D') and G C°°(T). Testing the equation of 

d 

the eigenfunction (f){k -f ijds^x) with -0 and , we derive 

oki 

dip 



{L{k ijdg) — Xj{k -1- i(3s))4>{k + i(3s, x)'ilj{x)^:—{k)dxdk = 0. 


D’ JT 


dh 


Observe that L(k -|- ijds)* = L{k — ijds) and 
integrate by parts to derive 


/ dL{k - i(3s) 

V dki 


dL{k -I- ifdg) 


dki 


We 


dip 


0=1 {{L{k + i/3s) - Xj{k + i/3s))(j){k + i/3s,x),i!{x))^2(j)-^{k)dk 

dip 


' D' 


(93) 


' D' 


ID' 


’ D' 


[4>{k i/3s, x), [L{k - i(3s) - Xj{k + i/3s)) ip{x)) 


LAD dki 


{k)dk 


d(j) 

dki 


{k + i/3s, x), {L{k - i(3s) - Xj{k + i/3s)) 'ip{x] 


L2(T) 


ip{k)dk 


j.n , ^ dL{k-if3s) ^ 

(p{k + ifds, x), --- f:(x) 


dki 


dXj 


{k -f il3s)'d{x) 


ip{k)dk 


' D' 


d(j) 


{L{k + i/Ss) - Xj{k + i(3s)) -^{k F il3s,x),'ilj{x) 

OKi 


L2(T) 

ip{k)dk 


L2(T) 


< D' 


/ dL{k + i(3s) dXj 

' VV dki dki 


{k -1- i(3s) (j){k + ifds, x),'d{x] 


ip{k)dk. 


L2(T) 


(t>i{.k,x) 

G{k) 

H{k,x) 


d(j) 

dki 


{k + i/3s,x), 


L{k + i(3s) - Xj{k + i(3s), 
fdL{k-^i(3s) dXj^j ^ 
V dki dki ^ 


(j){k -\- i/3s, x). 


We introduce 
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By invoking the previous lemma, the Lipschitz continuity of the C”^’''^(T)-valued 
function 0(/c + i(3s^ ■) implies that the mapping k i—)■ H{k^ ■) must be Lipschitz as a 
C'™'(T)-valued function on D' for any m > 0. On the other hand, the iL^(T)-valued 
function •) is also Lipschitz on D'due to (P4). Hence, both •), 

and {H{k, •),-^)i 2 (T) are continuous on D' for any test function The continuity let 
us conclude from (195]) that for every k G D', •) is a weak solution of the equation 

(94) G{k)(j)i{k,x) = H{k,x). 

We interpret (IMl) in the classical sense since all the coefficients of this equation are 
smooth. Consider any ki,k 2 in D' and subtract the equation corresponding to ki 
from the one corresponding to k 2 to obtain the equation for the oscillation function 

4>i{kir) -Mhr)- 

G{ki){(j)i{ki,x) - (j)i{k 2 ,x)) = {G{k 2 ) - G{ki))(j)i{k 2 ,x) + {H{ki,x) - H{k 2 ,x)). 

Note that due to regularities of Xj, H and the fact that the differential operator G{k) 
depends analytically on k, we get 

\\H{ki, ■) — H{k2 , ■)\\h^{t) + \\{G{ki) — G{k2))4>i{k2, •)||j7"*(t) = 0{\ki — k2\), Vm G M. 

Combining this with the uniform boundedness in k of the supremum norms of all 
coefficients of the differential operator G(fci), we obtain 

•) — •)lli7"*(T) = 0 {\ki — ^ 2 !), 

by using energy estimates as in the proof of Lemma [9~41 An application of the Sobolev 
embedding theorem shows that d^(j)i{k, x) is continuous on D' x T for any multi-index 

/S- 

To deduce continuity of higher derivatives d^d^4’{k + i(3s ) (|a| > 1, \(3\ > 0), we 
induct on I a I and repeat the arguments of the |a| = 1 case. 

Finally, the last statement of this proposition also follows since all of our estimates 
hold uniformly in s. □ 

Observation 9.7. 1. The property (P4) is crucial in order to bootstrap regularities 
of eigenfunctions 0(A; ± i(3s)- 

2. If one just requires ± ifis) G G^{D x for certain m > 0 then the 
smoothness on coefficients of L could be relaxed significantly (see [T2ll23] L 

10. Concluding remarks 

(1) The condition that the potentials A, V are infinitely differentiable is an 
overkill. The Fredholm property of the corresponding Floquet operators is 
essential, which can be obtained under much weaker assumptions. 

(2) The main result of this article assumes the central symmetry (evenness) of 
the relevant branch of the dispersion curve Xj{k), which does not hold for 
instance for operators with periodic magnetic potentials [101128]. Note that 
the result of [22] at the spectral edge does not require such a symmetry. It 
seems that in the inside-the-gap situation one also should not need such a 
symmetry. However, the authors have not been able to do so, and thus were 
limited to the case of high symmetry points of the Brillouin zone. 
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(3) In the case when A is below the whole spectrum, the result of this paper 
implies the Theorem 1.1 in [21] for self-adjoint operators. 
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